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Abstract: This survey is a collection of various results and formulas by 
different authors on the areas (integrals) of five related processes, viz. Brow- 
nian motion, bridge, excursion, meander and double meander; for the Brow- 
nian motion and bridge, which take both positive and negative values, we 
consider both the integral of the absolute value and the integral of the 
positive (or negative) part. This gives us seven related positive random 
variables, for which we study, in particular, formulas for moments and 
Laplace transforms; we also give (in many cases) series representations and 
asymptotics for density functions and distribution functions. We further 
study Wright's constants arising in the asymptotic enumeration of con- 
nected graphs; these are known to be closely connected to the moments of 
the Brownian excursion area. 

The main purpose is to compare the results for these seven Brownian 
areas by stating the results in parallel forms; thus emphasizing both the 
similarities and the differences. A recurring theme is the Airy function which 
appears in slightly different ways in formulas for all seven random variables. 
We further want to give explicit relations between the many different similar 
notations and definitions that have been used by various authors. There 
are also some new results, mainly to fill in gaps left in the literature. Some 
short proofs are given, but most proofs are omitted and the reader is instead 
referred to the original sources. 
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1. Introduction 

This survey started as an attempt to better understand two different results, 
viz. the asymptotic enumeration of connected graphs with a given number of 
independent cycles by Wright [H; [5§| and the formulas for moments of the 
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Brownian excursion area by Louchard [3j|, Takacs [4j| and others, and the 
surprising connection between these two seemingly unrelated results found (and 
explained) much later by Spencer [i^l, who showed that the same sequence of 
constants appears in both results (see ([36]) below). The literature may seem 
confusing (to me at least), however, because different authors, and sometimes 
the same author in different papers, have used not only different notations for 
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the same constants but also several different related constants. I tlierefore set 
out to collect the various definitions and notations, and to list the relations 
between them explicitly. I was further inspired by the paper by Flajolet and 
Louchard [17[, studying further related properties of the distribution. 

After doing this, I realized that many of the results for Brownian excur- 
sion area have close parallels for the integrals of (the absolute value of) other 
related processes, viz. Brownian bridge, Brownian motion, Brownian meander 
and Brownian double meander, and also for the positive parts of Brownian 
motion and Brownian bridge (see Sections [20l - [25)) . In particular, there are sim- 
ilar recursion formulas for moments and similar expressions for double Laplace 
transforms involving the Airy function, but the details differ between the seven 
different Brownian areas. I therefore decided to collect various results for these 
seven different areas of Brownian processes together so that both the similari- 
ties and the differences would become clear. One important source of results and 
inspiration for this is the series of papers by Takacs 44; 45; 3; 47; 48; 46 
another, treating different cases in a unified way, is Perman and Wellner [36 
the course of doing this, I have also added some minor new results motivated by 
parallel results. For example, much of Section and the related Appendix [Bl 
which both generalize results by Flajolet and Louchard 17|, seem to be new, as 
do many of the explicit asymptotic results as a; — > 0. 

To faciliate comparisons and further studies in the literature, I have often 
listed the notations used by other authors to provide a kind of dictionary, and I 
have tried to give as many explicit relations as possible, at the cost of sometimes 
perhaps being boring; I hope that the reader that gets bored by these details 
can skip them. I have found draft versions of this survey very useful myself, and 
I hope that others will find it useful too. 

There are two main methods to prove the results on moments for the various 
processes given here. The first studies Brownian motion using some form of the 
Feynman-Kac formula, see for example Kac 2^ 3^; 31 1 , Louchard 34 1 , Takacs 
|49j and Perman and Wellner [ssj : this leads to various expressions involving 
the Airy function. See also Section [27| for a probabilistic explanation of the two 
slightly different forms of formulas that arise for the different cases. A related 
method is the path integral method used in mathematical physics and applied 
to Brownian areas by Majumdar and Comtet [35| . The second main method 
is combinatorial and studies the corresponding discrete structures, combined 
with singularity analysis for the resulting genera ting functions; see for example 
Takacs 0; El Ei IS; S; [s^ ; and Nguyen The [Ig. 

We have omitted most of the proofs, but we include or sketch some. In par- 
ticular, we give for completeness in Appendix [C] proofs by the Feynman-Kac 
formula of the basic analytic identities (a double Laplace transform that is ex- 
pressed using the Airy function) for all types of areas of Brownian processes 
studied here. 

There are, of course, many related results not covered here. For example, 
we will not mention integrals of the supremum process or of local time [26|, 
or processes conditioned on their local time as in [1]. See also the impressive 
collection of formulas for Brownian functionals by Borodin and Salminen 
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The Brownian excursion area is introduced in Section [2] and some equivalent 
descriptions are given in Section[31 The results on graph enumeration are given in 
Section[4l the connection with the Brownian excursion area is given in Section[5l 
and various aspects of these results and other related results on the Brownian 
excursion area are discussed in Sections IBHTSl Integrals of higher powers of 
a Brownian excursion are discussed briefly in Section 1191 The other types of 
Brownian processes are studied in Sections [^DH^ these sections are intended 
to be parallel to each other as much as possible. Furthermore, in Sections [M] 
and I25i results are given both for the positive parts of a Brownian bridge and 
motion, respectively, and for the joint distribution of the positive and negative 
parts. Sections [26] and [27] give more information on the relation between the 
seven different Brownian areas (the reader may prefer to read these sections 
before reading about the individual variables) , and Section [^5] compares them 
numerically. Finally, Section [29] discusses negative moments of the Brownian 
areas. 

The appendices contain some general results used in the main part of the 
text. 

The first few moments of the seven different Brownian areas studied here 
are given explicitly in [11] tables throughout the paper together with the first 
elements in various related sequences of constants. Most of these values, and 
some other expressions in this paper, have been calculated using Maple. 

We let B{t) be a standard Brownian motion on the [0,oo) with B{Q) = 
(usually we consider only the interval [0,1]). We remind the reader that a 
Brownian bridge B\^^{t) can be defined as B{t) on the interval [0, 1] conditioned 
on -B(l) = 0; a (normalized) Brownian meander B^^^t) as B{t) on the interval 
[0, 1] conditioned on B{t) > 0, t G [0, 1], and a (normalized) Brownian excursion 
Bex{t) as B{t) on the interval [0, 1] conditioned on both B{1) = and B{t) > 0, 
t G [0,1]. These definitions have to be interpreted with some care, since we 
condition on events of probability 0, but it is well-known that they can be 
rigorously defined as suitable limits; see for example Durrett, Iglehart and Miller 



14| . There are also many other constructions of and relations between these 



processes. For example, see Revuz and Yor 39, Chapter XII], if T > is any 
fixed time (for example T = 1), and gr '■= max{i < T : B{t) = 0} and 
(It ■= min{<: > T : Bit) — 0} are the zeros of B nearest before and after T, then 
conditioned on and dr, the restrictions of B to the intervals [0,(7t], [qt^T] 
and [gxTdx] are, respectively, a Brownian bridge, a Brownian meander and a 
Brownian excursion on these intervals; the normalized processes on [0, 1] studied 

— 1/2 

in this paper can be obtained by standard Brownian rescaling as grp Bligx), 
(T - gT)-^'^B{gT + t{T - gr)) and [dr - gT)-^'^B{gT + t{dT - gr))- (See 
also Section 1271 where we consider Brownian motion, bridges, meanders and 
excursions on other intervals than [0,1], and obtain important result by this 
result and other rescalings.) Another well-known construction of the Brownian 
bridge is i?br(i) = B{t) — tB{l), t G [0, 1]. Further, the (normalized) Brownian 
excusion is a Bessel(3) bridge, i.e., the absolute value of a three-dimensional 
Brownian bridge, see Revuz and Yor [s^. 
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While the Brownian motion, bridge, meander and excursion have been ex- 
tensively studied, the final process considered here, viz. the Brownian double 
meander, has not been studied much; Majumdar and Comtet 35| is the main 
exception. In fact, the name has been invented for this paper since we have 
not seen any given to it previously. We define the Brownian double meander by 
Bdmit) := B{t) — mino<„<i this is a non-negative continuous stochastic 

process on [0, 1] that a.s. is at a unique point r S [0, 1], and it can be regarded 
as two Brownian meanders on the intervals [0, r] and [r, 1] joined back to back 
(with the first one reversed), see Section [23l 



2. Brownian excursion area and its moments 

Let -Bex denote a (normalized) Brownian excursion and 

Jo 

the Brownian excursion area. Two variants of this are 



(1) 



A 2-V2Sex (2) 

studied by Flajolet, Poblete and Viola [l^ and Flajolet and Louchard [l3| (their 
definition is actually by the moments in ([5]) below), and 



(3) 



used in [21|;|23|;|22]. (Louchard ISJ uses ^ for our Bex; Takacs [4J; |45|; |4a |47|; |4 



uses a;+; Perman and Wellner [38| use A^^cur] Spencer [43| uses L. Flajolet and 
Louchard use simply B for Sox, but since we in later sections will consider 
areas of other related processes too, we prefer to use Box throughout the paper 
for consistency.) Flajolet and Louchard [17[ call the distribution of A the Airy 
distribution^ because of the relations with the Airy function given later. (But 
note that the other Brownian areas in this paper also have similar connections 
to the Airy function.) 

This survey centres upon formulas for the moments of Box and various ana- 
logues of them. The first formula for the moments of Bex was given by Louchard 
(34 | (using formulas in fss'l), who showed (using /?„ for IE Box) 



where 7^ satisfies 



7r = 



12r r(3r + l/2) 



r-l 



6r - 1 r(r + 1/2) 



^ i r\ r(3j + 1/2) 

^ r(j -I- 1/2) 



7r- j J 



r > 1. 



Takacs O El El [i^; SI gives the formula (using Mk for EB^^) 



(4) 



(5) 



(6) 
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where Kq = —1/2 and 

3fc — 4 

Kk = -^Kk-i + KjKk-j, k>l. (7) 

(See also Nguyen The , where E B^^ is denoted and equals 2^/'^-''/'^Kr.) 

We win later, in Section [1^1 see that both the hnear recurrence (H])-® and 
the quadratic recurrence (H])-© follow from the same asymptotic series in- 
volving Airy functions. We will in the sequel see several variations and analogues 
for other Brownian areas of both these recursions. 

Flajolet, Poblete and Viola 18] and Flajolet and Louchard Jjj give the 
formula 

(and use it as a definition of the distribution of A; they further use /x^*^) and 
/Xfc, respectively, for E^*^). Here flk is defined by Flajolet and Louchard [l7| by 
Hq := —1 and the recursion 

k-l /,\ 

= {3k - 4)knk-i +Y.{ ■) ^i^k-j, k>l. (9) 

In particular, fli = 1/2. 

It is easily seen that ([5])-® are equivalent to (H))-© and the relation 

Flajolet, Poblete and Viola [3| define further 

ujk := ^ = 2^+'Kk, fc > 0, (11) 
ljI 2^''-^iJk = —^^k = 2^''Kk, k>0. (12) 
By (fT2l) . the recursion ([9| translates to 

u;l^2{3k-i)u;l_, + Y,^H~j, k > 1; (13) 

with uJq := —1/2 and uil = 1. Note that oj^ thus is an integer for fc > 1. The 
numbers are the same as in Janson j2l| . (This is a special case of uj^i 



in [21|, but we will only need the case I = 0.) The sequence (w^) is called the 
Wright-Louchard-Takacs sequence in [ist - 

By a special case of Janson 2l|, Theorem 3.3], 

92-5fc/2 
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which is equivalent to ([6]) by ([3|) and (fT2|l . 
Further, ^ is equivalent to ([5]) and 

7fc = fc > 0. (15) 

It is easily seen from ([9]) that 2'''0fc is an integer for all fc > 0, and thus by ((TSl) . 
7fe is an integer for fc > 0. 



3. Other Brownian representations 

Let _Bbr denote a Brownian bridge. By Vervaat ^5^, -Bbr(') ~ miii[o.i] -Bbr has 
the same distribution as a random translation of -Box, regarding [0, 1] as a circle. 
(That is, the translation by u is defined as i?cx(L' ~ ""J)-) particular, these 
random functions have the same area: 



6cx:= / Sex(t)di= / B,At)At~ min B^,{t). 
Jo Jo 0<t<l 



(16) 



This equivalence was noted by Takacs j4J; |45| , who considered the functional 
of a Brownian bridge in ()16|) in connection with a problem in railway traffic. 



Darling l^] considered the maximum (denoted G by him) of the process 

Y{t) := Bbr(i) - / Sbr(s) ds, 0<t<l. (17) 
Jo 

Since -Bbr is symmetric, Sbr = ^^br, we obtain from ()16|) also 

6ex= maxBbr(i)- / Bbr(i)dt== max r(i). (18) 
o<t<i Jo o<t<i 

Thus Darling's G := maxj Y{t) [l^ equals Bex (in distribution). Y{t) is clearly 
a continuous Gaussian process with mean 0, and a straightforward calculation 
yields 

Cov(y(s),r(t)) = i(|t-s|-i)2-^, s,te [0,1]. (i9) 

It follows that if we regard [0, 1] as a circle, or if we extend F periodically to M, 
y is a stationary Gaussian process, as observed by Watson [56[. Furthermore, 
by construction, the integral J^^ Y{t) dt vanishes identically. 

4. Graph enumeration 

Let C(n, q) be the number of connected graphs with n given (labelled) vertices 
and q edges. Recall Cayley's formula C{n,n~ 1) — n"^^ for every n > 1. Wright 
(sil l proved that for any fixed fc > —1, we have the analoguous asymptotic 
formula 

C{n, n + k)r^ pfcn"+(3'=-i)/2 as n ^ oo, (20) 
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for some constants pk given by 

2(l-3fe)/2^1/2 

r(3fc/2 + i) "^' ^''^ 

with other constants ak given by ct_i = —1/2, co = 1/4, ci = 5/16, and the 
quadratic recursion relation 

3(fc + l) 

(Tfc+i = (Tfe + o-jCTk-j, k>l. (22) 

Note the equivalent recursion formula 

3A: + 2 A , , 

o-fe+i = — - — CTfe + f^jf^k-j, k > -1. (23) 

i=o 

Wright gives in the later paper [s^ the same result in the form 

(although he now uses the notation fk = pk] we have further corrected a typo 
in [5^ Theorem 2]), where di = 5/36 and 

fc>i. (25) 

See also Bender, Canficld and McKay Corollaries 1 and 2], which gives the 
result using the same dk and further numbers Wk defined by wq = 7r/V6 and 

(8/3)i/27r(fc- 1)! /27A:\'=/2 , , , 

r(3fc/2) (^) 

so that 

31/2 . e \'=/2 

^^ = 2^(m) ^^0- (2^) 

(Wright [11 and Bender et al. [i| further consider extensions to the case /c — > 00, 
which does not interest us here.) 
In the form 

2-(5fe+l)/23fe+l^l/2^| 

= r(3fc/2 + i) 



(|24l) holds for all fc > 0, with da = 1/6. 

Wright's two versions ^ [58] and ([25]) [Hi] are equivalent and 

we have the relation 

/3\ '=+1 

afc = (-j fcldfe, fc>0. (29) 
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Next, define Cfc, fc > 1, as in Janson, Knuth, Luczak and Pittel ^24j, §8]; Ck is 
the coefficient for tlie leading term in an expansion of the generating function 
for connected graphs (or multigraphs) with n vertices and n + k edges. (Note 
that Cfe = Cfco = Cfco in [2I, §8]. ct is denoted Ck,-3k in Wright [HI] and bk in 
Wright [5^.) We have by Wright (sl, §5], or by comparing (|25|) and (f32|) below, 

cfc - (^)'(fc-l)!dfc, fc>l. (30) 
From [24, (8.12)] (which is equivalent to Wright 58, (7)]) follows the recursion 

3rcr = i(3r - l)(3r - 3)cr-i + f ^ j(r - 1 - j)cjCr_i_j, r > 1, (31) 
where jcj is interpreted as 1/6 when j = 0. In other words, ci = 5/24 and 

r-2 

rc, = ir(3r-3)c,_i + |^j(r-l-j)cjc,_i_j, r > 2. (32) 
By ([291), dSni) is equivalent to 

3 

o-fc = -^kck, fc > 1 (33) 

(and also for fc = with the interpretation Ocq = 1/6 again). 
By 0, §§3 and 8], ^ holds for fc > 1 with 

2(l-3fc)/2^ 

^^•^ r(3fc/2) 



which clearly is equivalent to and by ([33]) and ([301 
Finally, we note that ([20| can be written 



C(n,n + fc-l)~pfe_in"+3'=/2-2 ^ > q_ (35) 
5. The connection 

The connection between Brownian excursion area and graph enumeration was 
found by Spencer [43] , who gave a new proof of ([55]) , and thus (PH)) , that further 
shows 

Pk-i = fc > 0. (36) 

See also Aldous 0, §6]. 

By © and (I^T]) we thus have 

afc_i = 2^Kk, fc > 0, (37) 

and it is easily seen that is equivalent to ([7]) . Similarly, Q is equivalent to 
(|23)) and 

f2fc = 2fc!CTfc_i, fc>0, (38) 
and dH) is equivalent to jH]) by Jl]), dSB]) and 
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6. Further relations 

Further relations are immediately obtained by combining the ones above; we 
give some examples. 

By dMl) and ini), or by comparing ^ and 

Lol^2"'ak-i, k>Q, (39) 

and thus also, by (fT2|l . 

ujk = 2o-fc_i, fc > 0. (40) 
By (l38l) and p3|) . or comparing ((3T|) and ([91), we find 

r^fc = 3(fc- l)/c!cfe_i, fc>2. (41) 

By (Hg) and gl|), or by ^ and ([Ml), or by ^ and (US]), 



^^ = 22'=-i3(fc-l)cfe_i, k>2, (42) 
tjfc = 3(/c- l)cfe_i, A:>2. (43) 



By ^ and (p?!). 



/ 4 \ K, 



By (go]) and (l29|), or by ^ and (gi 

'3\fc 



u;, = 2(-) (A:-!)! 4-1, k>l. (45) 
By dSI]) and dM]), 

22-3fc/2^ ^ 22-7fc/2y^ 

r((3fc - 1)72)'^'=^' ~ r((3A: - l)/2)" 



Note further that can be written 



22-3fc/2 

^^-^ = r(3(fc-i)/2) ^^-^' '^2- 

By (Ull) and gH), or by ©, dM]) and gT]), 
as claimed in Janson [i^, Remark 2.5]. 



S. J anson/ Brownian areas and Wright's constants 



90 



7. Asymptotics 



Wright 'S!?! proved that the Hmit \m\k~,ryo dk exists, and gave the approximation 

0. 159155. The hmit was later identified by Bagaev and Dmitriev [Jl as l/(27r), 

1. e. 



2n 



as fc — > cxD. 



See [24|, p. 262] for further history and references. 
It follows by Stirling's formula that for Wk in ([^S)) . 



Hence, by ^ 
and, equivalently, 



Pk 



Wk ^ I 
31/2 



as fc — > 00. 



k/2 



pk-1 



27ri/2(i2fc) 



as fc — > cx). 



k/2 



\12kJ 



as k 



(49) 

(50) 
(51) 

(52) 



By (p6l) follows further, as stated in Takacs 44; 4^ 3 47; 4^, see also Janson 
and Louchard [25j, 



k \fe/2 



as fc ^ 00, 



(53) 



and equivalently, as stated by Flajolet and Louchard [17| and Chassaing and 
Janson respectively. 



V3e 



as fc ^ 00, 

E^fc _ 2i/23fcf^V as A; ^00. 
\3e/ 



3eJ 
k \fe/2 



By daO]) and (091), or by 0, Theorem 8.2], 

1 / 3\ 
^^-2^(2) 

Further, by ^ and (05]), 

1 /3\'=+i 
^^^■^2^12; 

and by (gll) and ([49]), as stated in Takacs [ii: lisi: liel: li?!: lij 

1 /3\fc, 



as /c ^ 00. 



as fc ^ cx), 



Kk 



27r V4 



(fc-1) 



as fc 



(54) 
(55) 

(56) 

(57) 

(58) 
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Jli^cx — 


1 
i 


n? R 

Ic-Ocx — 


V2^ 
4 


F k2 _ 5 

i5„^ — 

12 


l^'^ex — 


15V2^ 
128 


EOex — 


221 
1008 


E^" = 


1 


EA = 




1 n 

3 


= 


4 


Ey^* = 


"es" 


70 — 


-1 


71 = 


9 


72 = 405 


73 — 


65610 


74 = 


21749715 


Ko = 


1 

2 


i^i = 


1 

8 


X2 = A 

o4 


Ks = 


15 
128 


K4 = 


1105 
4096 


Oo = 


-1 


ni = 


1 

Z 


4 


^3 = 


45 

A 

4 


C4 = 


3315 

1 ft 

Id 


UJQ = 


-1 


ui = 


1 

o 
Z 


5 

o;2 = - 

o 


1^3 = 


15 

Q 
O 


0J4 = 


1105 

iZo 


UJq = 


1 

~2 


^1 = 


1 


= 5 


= 


60 


1^1 = 


1105 




1 


po = 




5 

''^ = ii 


P2 = 


5^2^ 


P3 = 


221 


p-1 = 


4 


256 


24192 


(T_l = 


1 

~2 


(70 = 


1 

4 


5 

"^ = 16 


cr2 = 


15 
16 


0-3 = 


1105 
256 






do = 


1 
6 


36 


d2 = 


5 

36 


da = 


1105 
7776 










5 

~ 24 


C2 = 


5 

16 


C3 = 


1105 
1152 










Table 1 











Some numerical values for the Brownian excursion and for graph enumeration. 

Similarly by dS?]) and (011), dMl) and 
1 



— 6%k~l)l as A: ^00, (59) 
i(^)V-l)! asfc^cx), (60) 



8. Numerical values 



Numerical values for small k are given in Table [TJ See further Louchard 34 1 

24, 



(ES^^), Takacs [4J; |47|; |48l [EB^^, Kk), Janson, Knuth, Luczak and Pittel 
p. 259 or 262] (cfe), Flajolet, Poblete and Viola [18i, Table 1 and p. 503] (E^*^, 
rife, ujk, ^*k), Flajolet and Louchard [13, Table 1] (E^'', ilfc), and Janson [2l|, p. 
343] {ul). 
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9. Power series 

Define the formal power series 

oo 

C{z):^^crz\ (62) 



By fi^, §8 and (7.2)], we then have 



Y,erz\ (63) 



r=0 



where 

(6r)! _ /3\'^r(r + 5/6)r(r+ 1/6) _ /3\'^(5/6)"(l/6)^ 



(i) 



25'^32'-(3r)! (2r)! \2J 27rr! 
^ r(3r + l/2) (64) 
r(r + l/2)r!' 

(The last formula follows by the triplication formula for the Gamma function, 
or by induction.) We have Cq = 1, ei = 5/24, ez = 385/1152, eg = 85085/82944. 
The constants are the coefficient for the leading term in an expansion of the 
generating function for all graphs (or multigraphs) with n vertices and n + k 
edges, see [24, §7] (where is denoted 6^0)- 



10. Linear recursions 

From ([63]) follows the linear recursion, see [13, §8], 

1 ''"^ 

Cr=er ^jCjCr^j, r > 1, (65) 



where are given explicitly by ([M)) . 

By (HT]) together with ([M)) and simple calculations, is equivalent to 

_ 12r r(3r+l/2) /r\ r(3j + 1/2) 

''""-^ r(. + i/2)-^^U r(, + i/2) i« ^^^^'(6^) 

given by Flajolet and Louchard [17], and, equivalently, see (fT5|l . 
_ 12r r(3r+l/2) /A r(3j + 1/2) 

- 6r - 1 r(r + 1/2) ^7 r(j + 1/2) " ^^^> 

given by Louchard [SJj; the latter was already given as ^ above. Changing the 
upper summation limit we can also write these as 

,0.0 _ 6. + 1 r(3r + 1/2) - /A m + 1/2) 

'''''■"6r-l r(r + l/2) A. ^ r(j + 1/2) ^ 



S. J anson/ Brownian areas and Wright's constants 93 

6^ r(3r+l/2) - /A r(3j + l/2) 

6r - 1 r(r + 1/2) j^^ \j) V{j + 1/2) ^^"^ ' " ^''^ 



By (fTUl) , these are further equivalent to the hnear recursion in Takacs 4J; |45 1 

6r + 1 

" 2(6r- 1) ""^ " 51 "j^'-^-J' ^ 1' C^O) 



where, cf. Ifil 



r(3j + 1/2) 
r0- + i/2)j! 



^-36-^ ;^^-;^;^, = 2-^e„ ,>o. (n) 



(Perman and Wellner [38] and Majumdar and Comtet [35| use jj for aj.) 
11. Semi-invariants 



Takacs [44| considers also the semi-invariants (cumulants) A„ of Bex, given re- 
cursively from the moments Mk '■= E^S^^ by the general formula 

A„ = Mn-^r 1 A4A„_fc, 71 > 1. (72) 

k=l ^ ^ 

For example, Ai = EiJex = ^/ttJS, A2 = Var(Kcx) = (10 - 37r)/24 and A3 = 
EB^^ - 3EB2^EBox + 2(EBex)^ = (Stt - 25)\/2^/128. Since these do not have 
a nice form, they will not be considered further here. 



12. The Airy function 

There are many connections between the distribution of Bex and the transcen- 
dental Airy function; these connections have been noted in different contexts by 
several authors and, in particular, studied in detail by Flajolet and Louchard 



17[. We describe many of them in the following sections. 

The Airy function Ai(x) is defined by, for example, the conditionally conver- 
gent integral 

1 f°° 

Ai{x) := - / cos(i^/3 + xt) dt, -00 < x < 00, (73) 
Jo 

see e.g. Abramowitz and Stegun 0, §10.4] or Lebedev [H, §5.17], which also 
give the following representations using standard Bessel functions 

Ai(^) = ^(i)'^'^i/3(^): ^>0, (74) 

^1/2 , ,0^3/2, ,0 3/2,, 

Ai(x) = ^ (/_i/3 (^) - /1/3 (^) ) , X > 0, (75) 



S. J anson/ Brownian areas and Wright's constants 



94 



Ai(-.) = — (^-1/3 (^) + ^1/3 (^) j , - > 0. 



(76) 



Ai(a;) is, up to normalization, the unique solution of the differential equation 
Ai"(a;) = a;Ai(a;) that is bounded for x > 0. Ai extends to an entire function. 
All zeros of Ai(a:) lie on the negative real axis; they will appear several times 
below, and we denote them by Cj = — |aj|, j = 1, 2, . . . , with < |ai| < |a2|, . . . . 
In other words, 

Ai(a,)=Ai(-|a,|) = 0, J -1,2,... (77) 

We have the asymptotic formula aj ~ — (37r/2)^/'^j^/^, which can be refined to 
an asymptotic exp ansion 0, 10.4.94], [13]. (Abramowitz and Stegun 0, 10.4.94] 



and Louchard (34| use aj as we do, while Takacs [441 : l45t l46t 1471: l48l : ISOl ] uses aj for 
\aj\; Fla .iole t and Louchard [13] and Majumdar and Comtet jSSlj use aj for \aj\; 
Darling [l^] denotes our \aj\ by cr^; his aj are the zeros of Ji/-i{x) + J_ 1/3(2;), 
and are equal to iJaj]^''^ in our notation, see (|76|l.') 

We will later also need both the derivative Ai' and the integral of Ai. It seems 
that there is no standard notation for the latter, and we will use 

/•oo -1 

kl{x) := / Ai{t) di = - - / Ai(i) dt, (78) 
using M{x)dx = 1/3 10.4.82]. See further AppendixEl 
13. Laplace transform 

Let ipcx{t) Ee~*^™ be the Laplace transform of Bex', thus tpcxi—t) '■= Ee*^"'' 
is the moment generating function of Bex- It follows from (|53p that tpcx{t) exists 
for all complex i, and is thus an entire function, with 

V^ex(t)=EE^cx^- (79) 

Uhx{t) is denoted (j){t) by Louchard l33l. G (2~^/'^t) by Flajolet and Louchard 
|17| and '^eit) by Perman and Wellner |38|.) 

Darling [I4] found (in the context of (|18p) the formula, with \aj\ as in ([77|) . 

00 

i/,^^(t) = %/2^t^exp(-2-i/3|^^. 1^2/3^^ ^ > (80) 



also found by Louchard [3J] using the formula, for a: > 0, 

Jo ^ '^"""^ 'V2^~ Ui(2i/3^) Ai(0)j' ^''^ 
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proved by Louchard [s^]- We give a (related but somewhat different) proof of 
(|5T|) in Appendix IC.4I The relation ([50]) follows easily by Laplace inversion from 
([81]) and the partial fraction expansion, see 

Ai(z) Ai(0) j^\,^a,^aj- ^ ' 

A proof of ([80l) by methods from mathematical physics is given by Majumdar 
and Comtet [35l |. 

Taking the derivative with respect to x in ([5T|) . we find, for x > 0, 



2 



r e--^ a^/^)^ - -2V3^ (^ Ai'(2V3,) X _ / 3 Ai'(2V3.) \ _ 

(83) 

or, by the changes of variables x ^ 2 ^^^x, t ^ 2^/'^x, as in [34 

7o dx V Ai(x) y V Ai(a:) / 

14. Series expansions for distribution and density functions 

Let Fcx(a;) := P(;Bex < a^) an d .f ex(x) := -F'(a;) be the distribution and density 
functions of Be^. (Takacs [H; Hi Eg [il; Sll uses W and W for these, while Fla- 
jolet and Louchard use W and w for the distribution and density functions 
of ^ = 2^/^Bex; thus their Wix) = F^^{2-^/^x) and w{x) ^ 2-^/^ f^^{2-^/^x).) 
Since Bex > a.s., we consider in this section x > only. 

Let, as above, aj denote the zeros of the Airy function. Darling 121 found by 
Laplace inversion from (j80p (in our notation) 

oo 

Fex(x) - 2V^J2 |aj r'/'p(V2|a, 1-3/2^), (85) 
i=i 

where p is the density of the positive stable distribution with exponent 2/3 , 
normalized to have Laplace transform exp(— t^/^); in the notation of Feller [l5|, 
Section XVIL6] 

pix)^p{x;l-l) = -1 £(-1)^-1 sin(2W3)^^^^^§^x-''=/^ (86) 

fc=i 

By Feller [H, Lemma XVIL6.2] also 

pi:,)^:c-'/^p{a:-y^;l^l), (87) 

where p{x; |,-i) = sin(2fc7r/3) x^ is the density of a 

spectrally negative stable distribution with exponent 3/2. (Takacs defines 
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the function g{x) — 2^^''''p(2^^/'^x; |, — i); this is another sueh density function 
with a different normahzation. Takacs [49| uses g for p.) The stable density 
function p can also be expressed by a Whittaker function or a confluent 
hypergeometric function U 0, 13.1], [H, 9.10 and 9.13.11] (where U is denoted 



-3\i/2 

TT/ 

24/3 

33/27ri/2 



X ''/^expl 



V 27x2/ V6'3'27a;2y 



Hence (US]) can be written [13; SI; [H; IH El] 

y6v^ 2/3 4 \ 



with Vj = 2|aj|3/27a;2, which leads to [Ulil IH IH; 133; Ss 



2V6 



(88) 
(89) 

(90) 
(91) 



(92) 



15. Asymptotics of distribution and density functions 

Louchard [34] gave the two first terms in an asymptotic expansion of fcx{x) and 
the first term for Fc:^{x) as x — s- 0; it was observed by Flajolet and Louchard 



171 ] that full asymptotic expansions readily follow from (for f^^; the result 
for i^ex follows similarly from (j90p . or by integration); note that only the term 
with j — 1 in (|92D is significant as x —^ and use the asymptotic expansion of 
U g iven in e.g. [3^, (9.12.3)]. The first terms are (correcting typos in |3J] and 
[U), as a; 0, 

/e.(x) ^ e"^ (^jail^/^-^ - ilail^/^x-^ ~ iflail-^/lx-i + ...), (93) 

FUx) - (fjail^/^x-^ + i JaiJ-^/^ _ ^laiJ-^/^a;^ + ...). (94) 

As in all other asymptotic expansions in this paper, we do not claim here that 
there is a convergent infinite series on the right hand side; the notation (using 
~ instead of =) signifies only that if we truncate the sum after an arbitrary 
finite number of terms, the error is smaller order than the last term. (Hence, 
more precisely, the error is of the order of the first omitted term.) In fact, the 
series in ([M]) and ([M]) diverge for every a; > because the asymptotic series for 
U does so. 
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For X — > CX), it was observed by Csorgo, Shi and Yor ll| that ([53|) imphes 

- lnP(6ex > 2:) = -ln(l- Fex(a;)) 6a;^ a; ^ 00. (95) 

Another proof was given by Fill and Janson ^IG*]. Much more precise results 
were obtained by Janson and Louchard [25|| who found, as a; ^ 00, asymptotic 
expansions beginning with 

^2V6,^_e^2/ 1 2 5 4 25 g 



/cx(a;) T^x^e-'"' 1 - - x"-* a;"'' + . . . , (96) 

^ V 9 1296 46656 I ' v 7 

6V6 _6a;^ 1^1 1 ^-2 1 ^-4 7 g 



1-Fe.(.)^^.e- (1-36-^-^-648^ - 46656^ • ^''^ 



16. Airy function expansions 

Following Louchard [sJl , we expand the left hand side of ([5T|) in an asymptotic 
series as a; ^ 00, using ([75)) . and find after replacing 2^/'^a; by a: the asymptotic 
series 

~ _2;i/2 + g(_l)fcE6,\^^^fc^2fc/Vi-3'=)/2. (98) 

Note that the infinite sum in ([95)1 diverges by ([55)1 : cf. the discussion in Sec- 
tion [ini For the proof of , we thus cannot simply substitute (|7^ into ([5T|) ; 
instead we substitute a truncated version (a finite Taylor series) 

Vex(i) = E^^cx^ + 0(i^+^), (99) 
A:=0 

where N is finite but arbitrary. In similar situations for other Brownian areas 
in later sections, for example with (|133p below for the Brownian bridge, this 
yields directly a sum of Gamma integrals and an asymptotic expansion of the 
desired type. In the present case, however, we have to work a little more to avoid 
divergent integrals. One possibility is to substitute ([99]) into the differentiated 
version ([84| . evaluate the integrals, subtract x~^/'^ from both sides, and integrate 
with respect to x. Another possibility is to write the left hand side of (|8ip as 



e 



---1)^- r(v.cx(t^/^) 



V27rF Jo V27rt3 



+ / e--(^ex(t^/^)-l)^ (100) 



and subsitute ([99]) into the third integral, noting that the first integral is —\f2x^l'^ 
and the second is a constant (necessarily equal to 2^/'^Ai'(0)/Ai(0)). 
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Using jU, we can rewrite ([98]) as 

k=0 

This is also, with a change of variables, given in Takacs [i^. Flajolet and 
Louchard 17[ give this expansion with the coefficients (— f )'^rife/(2'^fc!), which 
is equivalent by (fTU]) . They give also related expansions involving Bessel and 
confluent hypergeometric functions with coefficients flk/k\. An equivalent ex- 
pansion with coeficients 3{k — l)ck^i, which equals flk/k\ by (dJ), was given by 
Voblyi [El, see (8.14) and (8.15)]. 

The Airy function and its derivative have, as a; — s- oo, the asymptotic expan- 
sions P, 10.4.59 and 10.4.61] (we write c'j, and dj. instead of Ck and dk to avoid 
confusion with our Cfe and dk above) 



Ai(x) ^ ^.-V4,-2.-/3^(_i).,. ,-3./2^ (102) 



OO 



Ai'(x) ~ _^xV4e-2.^'V3^(_i)ferf. (|)'=^^3fc/2^ (103) 

fc=0 



where 



r(3fc-|-i) 6k + l , , , , 

' dk = -^, 74, ^>0. (104) 



54fefc!r(fc-h i)' 6fc- 

These asymptotic expansions can also be written using or ak, since, by ([64|) . 
dni) and (fTUll . 

3'=4=efe, (105) 

(3/2)'=4=afe, (106) 
6A- + 1 

(3/2)'^4^4^:=-^^afc. (107) 

We thus obtain from (|10ip - (|103p , after the change of variables z = —x~'^l'^ /* 0, 
the equality for formal power series 



(108) 



fc=0 



These asymptotic expansions can also be written as hypergeometric series. 
We have, as is easily verified, the equalities for formal power series 



^c'fcz"=2i^o(5/6,l/6;z/2), (109) 

OO 

^4z^-2^^o(7/6,-l/6;z/2). (110) 



fe=0 
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Hence, (|102p and (|103p can be written, for x oo, 

Ai(.)^^.-V4,-..-/3^^„(5^1^_J_^^ (111) 

Ai'(.)^-^.VV--/3,Fo(I,-l;-^), (112) 
and ()108|) can be written 

>^ I. 2^0(7/6,-1/6; 3z/4) , , 

V 2KkZ^ = - '^J ' — f^. 113 

t^o 2Fo(5/6,l/6;3z/4) ^ ^ 

Flajolet and Louchard [l7| give the equivalent formula, see ((1 



f^^ _ 2^0(7/6,-1/6; 3W2) 

^' " 2Fo(5/6,l/6;3«;/2) ■ ^''^^ 

By the asymptotic expansion for Besscl functions [s^, (5.11.10)] 

h{x) ^ e''{2iTx)-^'^2Fo{ll2 + u,l/2~v]l/2x), x 00, (115) 
this can, as in Flajolet and Louchard [l3], also be written (with arbitrary signs) 
-^±2/3(3^/3) x^r, "^'^ V" ~k n1R^ 

The coefficients in this asymptotic series can be rewritten in various ways by 
the relations in Sections [2] and [6l for example by (|4T|) as 3(fc — l)ck-i for k >2, 
which yields ^ (8.15)]. 

If we in (jlOSp multiply by the denominator and identify coefficients, we obtain 
the linear recursion ([70)1 and the equivalent recursions ([S5)) - ([M)) : note that = 
3'^c; and ak = (3/2)*^ c'^ by (HHS]) and (fTTOl) . 

On the other hand, by differentiating (|10ip (which is allowed, e.g. because 
the asymptotic expansion holds in a sector in the complex plane) , 

-"'--™-9f-(f^)^ 

00 

^(-l)'=2iffea;i/2-3fc/2 ] . (117) 



fc=0 

/ ^ \ 2 



\fe=0 

This gives the equation for formal power series 



^(3A;-l)X,z'=+i = l-^2X,zM , (118) 



fc=0 \k=0 

which is equivalent to the quadratic recursion ([7|). 
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17. A continued fraction 

As noted in [23, (8.16)], (|116p yields (after substituting z = \/x) the asymptotic 
expansion 

■1 oo ^ oo 

where the continued fraction has denominators (6n + 2)z. 



18. Moments and Airy zeroes 



Define, following Flajolet and Louchard [17 1 



A(s) := ^ |aj|"", Res > 3/2. (120) 



Since \aj\ (37rj72)^/'^ P, 10.4.94] the sum converges and A is analytic for 
Res > 3/2. (Flajolet and Louchard [17| call A the root zeta function of the Airy 
function Ai.) 

By ([80l) . for Res > and with absolutely convergent sums and integrals, 



/•oo j*oo oo 

/ t"^Vex(i)dt = %/2^ / t''+i^exp(-2"i/3|aj|t' 

-'0 "'0 .1 



^ .=1^° (121) 



= 3Y^ gf^+l)/2 ^ |^^|-f3s+3)/2p|^ 3s + 3 j 

.3V..-r(4l)A(5£±l). 

In particular, this is finite for all s > 0, and thus Lemma IB. II in Appendix IB] 
applies and shows that Sox has negative moments of all orders. Since all (posi- 
tive) moments too are finite, ES^^ is an entire function of s, and by (jl2ip and 
P35[) . A extends to a meromorphic function in C and, as shown by Flajolet and 
Louchard |17j, 

^K.= '-^T(l^V{'-^\ ..C. (122) 



r(-.) 

Since the left hand side has no poles, while the Gamma factors are non-zero, 
A((3 — 3s)/2) can have poles only at the poles of r(— s), i.e. at s = 0, 1, ... ; 
hence, A-{z) can have poles only at z — |, 0, — |, —3, — |, .... Moreover, ^B^^ 
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is non-zero for real s; hence r((3 — 3s)/2)A((3 — 3s)/2) and r(— s) have exactly 
the same zeros (i.e., none) and the same poles. Consequently, A(z) has simple 
poles at z = |, — |, — |, . . . , but is finite and non-zero at 0, —3, —6, . . . (where 
r has a pole); furthermore, A is zero at the other poles of F, i.e. at —1, —2, 
—4, —5, —7, .... This was shown by Flajolet and Louchard [l7| (by a partly 
different approach, extending A(s) by explicit formulas). 

As shown in Flajolet and Louchard 17[, as a simple consequence of for 
\z\ < |ai|, 

Ai'(z) Ai'(0)^^^^^^^^^_^^,^ ^^23) 



Ai(z) Ai(0) 



fc>i 



and thus the values of A(s) at positive integers s = 2,3, .. . can be computed 
from the Taylor series of Ai, given for example in 10.4.2-5]. This and (|122p 
gives explicit formulas for the negative moments when s is an odd multiple 
of 1/3, including when s is an odd integer; see Flajolet and Louchard [17| . 
Alternatively, these formulas follow from p40p in Appendix |B] For example, 

9l3/6o2/3_5/2 

Ee-1/3 ^ ^ = 2-i^/63i9/6^"5/2p(2/3)5 ^ 1.184, (124) 

F(l/3)"' 

EB-^' = 3V2^ f 1 - ^!!^Lml) « 1.693, (125) 



47r3 

see further Section [29l We have here used the standard formula 

r(l/3)F(2/3) = — ^ = ^. (126) 
sm(7r/3) V3 

Note also that if —1 < Res < —1/2, then by ([5T|) . Fubini's theorem and 
((3351) in Appendix El 

21/3 r,.-i / Ai'(2V3,) Ai'(0)\ 
^ I ^ [ Ai(2i/3^) Ai(0) J 

'"":r^-^(e-"* - l)^cx(t'/')^^dx 



"'0 



V2^ 



r(.) 2 r„-2./3-./s^,_,„,.<i" 



V27r 3 7o w 
2r(s) / 2.s + l\ /9.,iw, 



3V27r 



By (HHl), this equals 



2i/3-s/3p(g)P(^ _ ^^^^ _ ^) ^ 2(i-^)/3 ^ A(i _ s), (128) 

sin(7rs) 



as shown directly by Flajolet and Louchard (17|. 
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19. Integrals of powers of i?ex 

Several related results are known for other functionals of a Brownian excursion 
or other variants of a Brownian motion. We describe some of them in this and 
the following sections, emphasizing the similarities with the results above, and 
in particular linear and quadratic recurrencies for moments and related formulas 
for generating functions. 

The results above for moments of Sex = /q Bcxit) dt have been generalized 



to joint moments of the integrals Bc^itY dt, £ = 1, 2, . . . , by Nguyen The 37| 
(mainly ^ = I and 2) and Richard |41|] (all ^ > 1). These papers show that ^ 
and d?]) (or equivalent formulas above) extend to these joint moments, with a 
quadratic recursion with multiple indices. We refer to these papers for details, 
and remark only that if we specialize to Kox {(■ = l)i then ^^.s in Nguyen The 



37[ specializes to 

Arfi^r'Kr^dr-i, r > 0; (129) 



similarly, in Richard 41|, 



/fc,o,...,o = 23-2fciffe, A: > 0, (130) 
mfc^o,....o = 23'^/2ESe^^ = Eyl^ k>Q, (131) 



so Xi^n — * -A, and ?/fe,o,....o = EB^^. The recursions in both 37] and 41 1 
specialize t o ([71 ) (or an equivalent quadratic recursion). 

Janson |2l| has given further similar results, including a similar quadratic 
recursion with double indices, for (joint) moments of Sex and another functional 



of a Brownian excursion; see also Chassaing and Janson 23|. 



20. Brownian bridge 

Let Sbr := /q^ \B\y^{t)\dt, the integral (or average) of the absolute value of a 
Brownian bridge. (Shepp [i^ and Ric ejioj use ^ for Sbi-; Johnson and Killecn 
[28] use L; Takacs uses a. Cifarelli [9| uses C{x) for the distribution function 
of 23/2Sbr.) Comtct, Desbois and Texier [13] give also the alternative represen- 
tation 

Sbri- / \BUt)\dt= / lBbr(i)-Bb,.(l-t)ldi, (132) 
Jo Jo 

which follows because Sbr(t) - Sbr(l - t) = -Bbr(2<), t e [0, 1/2]. 

Cifarelli [9] and (independently) Shepp ^4^] found a formula for the Laplace 
transform V'br(i) '■— Ee^*'''''. (See Appendix IC . 2 1 for a proof) An equivalent re- 
sult, in physical formulation, was proved by path integral methods by Altshuler, 
Aronov and Khmelnitsky see Comtet, Desbois and Texier [l3]. Shepp's ver- 
sion is (with Vbr denoted <f) by Shepp), cf (|5T|) : 

r e--V^br(V2i3/2)i^ = -V^^, X > 0, (133) 
Jo Vt Ai [xj 
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while Cifarelli's version, which actually is stated using the Bessel function iiTi/a, 
cf. ([7^ . can be written 

for arbitrary u > 0. Clearly, (|133p is the special case u = V2 of (|134p . and 
the two formulas are equivalent by a simple change of variables. The moments 
of ,Bbr the are obtained by asymptotic expansions, cf. Section [111 Cifarelli 0] 
considers (|134p as u ^ and Shepp [4§| considers (|133p as x — > oo; these are 
obviously equivalent. Introduce Z3„ defined by 



r((3n+l)/2)^ 



'■K.^^TT:;:r-TTT^^D^^ " > 0- (135) 



(Takacs [43] denotes E^B^J^. by M*; Perman and Wellner [si] use L„ for our !)„, 
^0 for Sbr and At„ for ES{J^; Nguyen The dH denotes EKJ^^ by and uses 
Q„ for our D„ and for 2-i/2-«/2^^^ -j ^j^g^^ obtain from cf. 
and the discussion after it, 



Ai(a;) 
Ai'(a;) 



Y,{~lfD,x-^^''-''\ (136) 



fe=0 

and thus by pU^ . pUg]) . pTO]) . pUT]) . cf. pUS]) . as formal power series, 

VP... Er.o(3/2)H^^ _ Er=o«.^'- .,,7. 

Cifarelli [§| uses 

mfe 23'=/2e6^V = E(2''^/^Sbr)''^ (138) 
and writes the result of the expansion as 

'"'- r((3f/i)/2) '-""^'- '^o- 

with, as formal power series (we use a'j^ and 6'^ for aj. and 6^, in 9]) 



5] C,z^ = where < = (-3)'=c;, and 6^, = (-3)'=^;,; 



(140) 



this is by (|138p equivalent to (|135p and (I137P together with 

Cfe = {-2fDk, k>0. (141) 

Cifarelli ^ gives, instead of a recursion relation, the solution to (|140p by the 
formula 

k k 

i-l^Ck = E(-l)'"'<-pl^lp = E K-p\\B\p, (142) 

p=0 p=0 
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where \B\p is given by the general determinant formula 



\B\, 



1 b[ 
1 







K-3 K-2 



1 b[ 



this is a general way (when 6q = 1) to write the solution to 



(143) 



(144) 



Shepp [431 gives, omitting intermediate steps, the result of expanding p33p 
as (using for Bhr) 



where e„ satisfies the linear recursion 

_r(37i+i) " _ fn\6k+ir{3k + ^) 



n > 0, 



kj 6k-l r(/s+ i) 



n > 0. 



(145) 



(146) 



(Shepp writes e„ but we use e„ to avoid confusion with in (l64|) .l 
We have, e.g. by (fUS]) . (fT35l) and (fT4T|) . 

e„ = 36"n! = (-18)"n! C„, n > 0. 



(147) 



Hence, (|146p can, as in Perman and Wellner ^] (where further, as said above, 
an is denoted by 7„), be rewritten as 



y 7^(^^Dn-^, n>0, (148) 



which also follows from (|137p . 

Takacs [J^l found, by different methods, a quadratic recurrence analogous to 
([U and (III), namely ()135p together with = 1 and 



3n - 2 1 

= T ^ 9 D,Dn-t, n>l. 



(149) 



This follows also by differentiating (|136p . which yields, as a; — > oo 
^(_l)fe+i^^!_Lli:)^.2;-3fc/2-3/2 ^ _i < Ai(x)Ai"(x) 



fe=0 



Ai'(a;)2 
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-1+ E(-i)'^^-= 



.-3/c/2 



\k=0 



and thus 



(150) 



3fc + 1 



kZ 



k=0 



\k=0 



Note further that (fTOTj) and (fT36ll show that 



E^^ 

fe=0 



fc=0 



EZoKkz'^ 



and thus 



Y^D,Kk-^ = 0, k>l, 



(151) 



(152) 



(153) 



i=0 



which can be regarded as a linear recursion for D„ given i4r„, or conversely. 
Explicitly, 



2^Ai^„-», n>l. 



(154) 



4=0 



It follows easily from the asymptotics ([55)1 that the term with i = in (|154p 
asymptotically dominates the sum of the others and thus 



IV. 



2 /3n\" 
7m V 4e / 



as n — > oo, 



(155) 



as given by Takacs 47[ (derived by him from (|149p without further details). 
Consequently, cf. ([55)) . by (|135p . or alternatively by Tolmatz 5l| or Janson and 
Louchard ^25] , 

, „ V n/2 

as n — > oo. (156) 



E 



The entire function Ai' has all its zeros on the negative real axis (just as 
Ai, see Section [T2|) . and we denote them by a'j = ~Wj\' J — 1j2,..., with 
< la'il < 141 < . . . ; thus, cf. ([TTD, 



Ai'(a')-Ai'(-|a'|) = 0, j-1,2,. 



(157) 



(Abramowitz and Stegunjl 



use a'j as we do; Takacs 14 



10.4.95], Rice ^Q,] and Johnson and Killeen 28 1 
49| uses a'j for our \a'j\; Flajolet and Louchard [17 1 
use aj for |a^|; Kac |29| and Johnson and Killeen [28[ use Sj for 2~-'^/^|a^|.) 

The meromorphic function Ai(z)/Ai'(z) has the residue Ai(a^)/Ai"(a^) = 
1 /a' at a'j , and there is a convergent partial fraction expansion 



oc ^ oo ^ 

.^-^ \n'\( 7 A- \n'\\ ' 



Ai(z^_ 

Ai'(z) ^a;.(z~a;.) 



(158) 
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Inversion of the Laplace transform in ()133p yields, as found by Rice [4o| . 

oo 

^brW =2-i/V/2^|a;.r'i'/'e-'"'''"^'*''', * > 0. (159) 
j=i 

Johnson and Killeen [2^ found by a second Laplace transform inversion the 
distribution function F\j,.{x) P(Sbr < x), here written as in Takacs 47| with 
u, := |a;.|V(27a;2), 

OO 

F^,{x) = 18-i/67ri/2a;-i^e-«^u7^/'Ai((3?/,/2)2/3) 

oo (160) 

^ 2-1/632/3^1/2^-1/3 y lg-|a;.|V(27x^)Ai(l8"2/3|„M2^-4/3y 



Numerical values are given by Johnson and Killeen [28| and Takacs [47 1. 

The density function /br = F{^^ can be obtained by termwise differentiation 
of l|160p . although no-one seems to have bothered to write it out explicitly; 
numerical values (obtained from (|159p ) are given by Rice [40|. 

Asymptotic expansions of fhr{x) and Fbr(a;) as x ^ follow from (|160p and 
P02p ; only the term with j = 1 in (|160p is significant and the first terms of the 
asymptotic expansions are, as a; — )■ 0, 



(161) 

Fb..(x) ^ (1 |al|-3/2 - II \a',\-^/^x' + iff lair^^/^x^ + O {x^)) . 

(162) 



For a; — > 00 we have by Tolmatz 5l| and Janson and Louchard [25| asymptotic 
expansions beginning with 



Nguyen The [37| extended some of these results to the joint Laplace trans- 
form and moments of Sbr and |i?brP dt. Specializing to Bhr only, his moment 
formula is 

with (correcting several typos in [s^]) -Bo,o = 1 and 

3r-2 ''"^ 
= ^^B,_i,o- i^B,,oS,_,,o, r>l. (166) 

i=l 
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4V 2 60 512V 2 '"'^ 720 

1 _ 14 21 ^ 76 
mn = 1 mi = —J-k = — 013 = — v/i" "^4 = — 

2 15 32 45 

eo = 1 ei = 9 62 = 567 £3 = 91854 64 = 28796229 

1 7 21 1463 

Do = 1 Di = - 1)2 = — -D3 = — ^^4 = 

4 32 64 2048 

Co = 1 Ci = — C2 = — C3 = C4 = 

2 8 8 128 

e J = 1 ej; = 1 §2 = 7 63 = 126 ej = 4389 

= 1 Dl =2 = 14 _DJ = 168 = 2926 

Table 2 

Some numerical values for the Brownian bridge. 



Clearly, these results are equivalent to p35p and (|149p . with 

Br,Q = TDr = {-lyCr: r>{). (167) 



Nguyen The (37| gives also the relation 

'A,flBr-^,0 = Q, r>l, (168) 



i=0 



which by (fT29)) and p67l) is equivalent to p33)) . Equivalently, 

Some numerical values are given in Tabled see further Cifarelli ^ {mn), 
Shepp [H (ES^^, e„), Takacs 47] (ES« , i?„), Nguyen The [H (ES;^' ^il^)- 
We define D* := 23"i:i„ and g; := 9-"e„ = 4"n! £>„ (considered by Shepp (H); 
note that these are integers. 

21. Brownian motion 

Let Bhm '■— Jq \B{t)\ dt, the integral of the absolute value of a Brownian motion 
on the unit interval, and let V'bm(i) := Ee^*^*"" be its Laplace transform. Kac 
[2^ gave the formula 

00 

V'bm(t) = t > 0, (170) 

where Sj = 2-'^/^\a'j \ and, with Kac's notation P(y) = 2i/3Ai(-2i/3y) , 
_ 1 + 3 J^' P{y) dy _ 1 + 3 L° Ai(2;) Ax 



i6,P{5,) 3|a'|Ai(a'; 



(171) 
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(Takacs [i^ uses Cj for Kj.) Using the function AI defined in ([78|) . (|17ip can 
be written 



AI(a;.) 
l«^|Ai(a^)- 



(172) 



In fact, see Takacs 4!| and Perman and Wellner 38], or (|362p in Appendix lC.il 
below, 

re"-'^b.n(%/2t3/2)di = -^, x>G. (173) 
Jo ' ' Ai (a;) 

The meromorphic function — AI(z)/Ai'(z) has residue — AI(aj)/Ai"(o^) — kj at 
a'j , and the partial fraction expansion 



AI(^) 

Ai'(z) ^z-a;. 



y^^y^^ ,eC. (174) 

r — n'. ^ r 4- \n'\ ^ ' 



An inversion of the Laplace transform in (|173p thus yields (|170p . 

Takacs [4§| found by a second Laplace transform inversion, this time of (|170p , 
the density function fhni{x) (by him denoted h{x)) of B^ra as 

oo 

hU^) - ^K,V2|a;.r3/2p(V2|a;r'/'2;) (175) 

-;^£».K)^"'-'''(s.^";)- c™) 

where p is the stable density in (f87|) . and w^- = 2|a^ p/27a::^ (denoted by Takacs 
Hi), cf. (ESI), dSIl), JMl). 

Asymptotic expansions of /bm(a^) and the distribution function i^bm(a;) as 
a: ^ follow from (|176p using the asymptotic expansion of U given in e.g. [13, 
(9.12.3)], cf. Section[T5l only the term with j = 1 in (|176p is significant and the 
first terms of the asymptotic expansions are, as x — > 0, 

f^^x) ^ ^e-^ (||a;i^/^x-^ + m-'/' - '^Kr/^x^ + ...), (177) 



F^^{x) ^ -^e-^ - '^Kr'^x' + '-^K^'^'x^ + ...). 

(178) 



For a; ^ 00 we have by Tolmatz [52| and Janson and Louchard [25| asymptotic 
expansions beginning with 

1 - Fb,„(i) ~ ^ x'^e-"''^ (l-Lx-' + ^x-' + ^..). (180) 
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Takacs j43] also found recursion formulas for the moments. Define L„ by 

O-n/2 \ 

(Takacs j49>] denotes ^B^^ by /i„; Perman and Wellner use for our 
L„, A for Sbm and for E6[,\„; Nguyen The uses for 2""/2^„.) An 
asymptotic expansion of the left hand side of (|173p yields, arguing as for (|98p 
and p36)) . 

AI(a:) g(_i)fci^^-3fe/2-i^ ^^^^ (182) 



Ai'(2;) 



fe=0 



Recall the asymptotic expansions (|102p and (|103p of Ai(a;) and Ai'(a;); there is 
a similar expansion of AI(a;) [4§|, see Appendix [XI 



oo 

,3/2 



fc=0 

where (3k, k > 0, are given by /Sq = 1 and the recursion relation (obtained from 
comparing a formal differentiation of (fTSS]) with (fTOS]) ) 

= + fc>l. (184) 

(Takacs fio^ further uses hk = (2/3)'=/3fc.) Hence, yields, using piI5)l and 

P07p . the equality (for formal power series) 

Multiplying by the denominator and identifying coefficients leads to the recur- 
sion by Takacs [i^, 

Ln = Pn + y^^ „. ^ CtjLn-j , n > 0. (186) 

Differentiation of (|182p yields, using Ai"(x) ~ a;Ai(x), 



Ai(x) xAi(a;)AI(a;) 
Ai'(a;) Ai'(a;)2 



g^_^^,+i3fc + 2^^^_3,/2_2^ (187) 



and thus by (|136p and (|182p the equality 

oo oo oo oo oi, I o °° OL 1 

- E ^'^-^ + E E = E ^i'.-'^-*-^ - E 

fc=0 fc=0 k=0 k=0 k=l 

(188) 
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= 1 




2 / 2 

sV^ 




3 

8 




/3o 


= 1 


/3i = 


41 
48 


/32 = 


9241 
4608 




Lo 


= 1 


ii = 


1 


L2 = 


9 

4 


is 


Lo 


= 1 


LI = 


8 


L'2 = 


144 


is 



630 V 71- 2560 
5075225 _ 5153008945 

663552 ~ 127401984 

263 ^ _ 2709 

"32" " ~ 64 

L;^ = 4208 

Table 3 

Some numerical values for the Brownian motion. 



LI = 173376 



which is equivalent to Lq = 1 a-nd the recursion formula 
3n — 1 

^ Ln-l = 2^ Dn-jLj, n > 1 



or 



i'n — ^ Ln^i - ^ Dn-jLj, 



n > 1. 



(189) 



(190) 



Takacs |49|] proves also the asymptotic relations, see further Tolmatz [52] and 
Janson and Louchard [2^ , 



1/2 



L„ ^ ^/3(-) . 



(191) 
(192) 



Some numerical values are given in Table [31 see further Takacs [49] and 
Perman and Wellner [sl]. We define L* := 2^"L„; these are integers by (I190p . 



22. Brownian meander 



Let B-mc '■— Jq \Bnie{t)\dt, the integral of a Brownian meander on the unit 
interval, and let V'me(i) := Ee"*^"" be its Laplace transform. 

Takacs [H^] and Perman and Wellner [s^ give formulas equivalent to 



-Vme(y2i3/2) 



dt Al{x) 



a: > 0; 



(193) 



see also Appendix IC.3I 
Define (5„ by 



\/^2- 



r((37i + 1)/2) 



n> 0. 



(194) 



(Takacs [5(1] denotes ^B^^ by M„ and V'me by Perman and Wellner 
Rn for our Qn, ^moan for Bmc and p„ for EyB^ie; Nguyen The [s^] denotes EB'^ 



use 
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by and uses for 2^/^^"/^Q„; Majumdar and Comtet [sl] denote ES,",, 
by a„ and Qn by Rn-) An asymptotic expansion of the left hand side of (|193p 
yields, arguing as for dM]), ([1^6)) and (fT82)) . 

^^f:(-l)'^g.x-3'^-/2-i/2^ (195) 

The asymptotic expansions (|102p and (|183p of Ai(a::) and AI(a;) yield, using 
(|106p . the equality (for formal power series) 

Differentiation of (|195p yields 

- 1 - ^^t]^^'^ - Y.i-ir'^-^Q,^-''"'-''^ (197) 
and thus, using (jlOip . the equality 

fe=0 k=0 fc=0 

which is equivalent to the recursion formula 

371 — 2 " 

2E^j<3n-j: «>1- (199) 

J=l 

with Qo = Ij proved by Takacs [50] by a different method, and the differential 
equation 

J C30 oo oo 

3z'— E Qfc^' + + 4 E ^fc^') E Qfc^"' = (200) 

fc=0 fe=0 fc=0 

also given by Takacs [s^ (in a slightly different form, with —z substituted 
for z). On the other hand, multiplying by J2T=o'^^'' CM]) and identifying 
coefficients leads to another recursion by Takacs SOg and Perman and Wellncr 

n 

Qn = Pn-J2 ^jQn^J' ^ (201) 

Furthermore, pI75|) . and (fTM)) yield the equality 

OO OO OO 

-2E^fe^'E^''^' = EQfe^^ (202) 

k=0 k=0 k=0 
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Qo 
Qo 



EBn 



- EBl^ = — 



59 
60 
59 
32 
118 



465 
512 
465 

"eT 

: 3720 



- EBt 



Table 4 

Some numerical values for the Brownian meander. 



5345 
3696 
80175 
2048 
160350 



and thus the relation 

n 

Qn = -2^KjLr,-j, n > 0. (203) 

3=0 

By (fT^ we obtain from pl^ also 

OO OO OO 

Y.DkZ^Y.Qkz" ^Y.^'^^"' (204) 

fc=0 fc=0 fc=0 

and thus the relation 

n 

L^ = Y,DjQn-j, n> 0. (205) 
Takacs [H^ proves the asymptotic relations, see also Janson and Louchard 

, / 77 \ n/2 

¥.BI,,^VZ^[-J , (206) 
Q„^V3(|-) . (207) 

Some numerical values are given in TablejH see further Takacs [H^. We define 
:= 23"Q„; these are integers by ((2031) (or by (fT99)) and (fT2|l ). 
The meromorphic function AI(z)/Ai(z) has residue rj := Al{aj) / Ai' (uj) at 
Oj, and the partial fraction expansion 

Al{z) ^ / I 1 \ 1 

where the first sum is conditionally convergent only and the second sum is abso- 
lutely convergent, since \aj\ x j^/^ and \rj \ x j^^/^. (Note that rj alternates in 
sign.) In fact, the second version follows from an absolutely convergent partial 
fraction expansion of AI(z)/(zAi(z)), using AI(0) = 1/3, and the first then fol- 
lows easily. A Laplace transform inversion (considering, for example, the Laplace 
transform of y^ipmciV^t^^^)) yields the formula by Takacs (50| 

OO 

V'„.e(0 = 2-W/3V^^r,e-2"'''l''^l*''', t> 0. (209) 
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A second Laplace inversion, see Takacs [50| . yields the distribution function 
Fmc{x) of Sine as, using Uj — |aj|^/(27x^) and Rj = \aj\rj = \aj\Al{aj)/ Ai' {a j), 



2/3\ 



Fn,o{x) = 18-i/67ri/2a;-i^i?^.e-".u7i/3Ai((3uj/2) 

oo (210) 

= 2-1/632/3^1/2^-1/3 ^^^.g-|a,|V(27.^)Ai(l8-2/3a2^-4/3)^ 
J = l 

Numerical values are given by Takacs [H^l- The density function /ino(a^) can be 
obtained by termwise differentiation. 

A physical proof of (|209p by path integral methods is given by Majumdar 
and Comtet [35| (denoting our rj by B(aj)), who then conversely derive (|193p 
from (|209)) . 

Asymptotic expansions of /mc {x) and i^mc (a;) as a; ^ follow from (|210|1 and 
(|102p ; only the term with j = 1 in (|210p is significant and the first terms of the 
asymptotic expansions are, as a; — > 0, 

Ue{x) - ne-^ (I \a,f/'x-' - A lail-i/^x-i + § laj-^/^ + O (x^)) , 

(211) 

F,ne(:r) ^ ne-^ (f |airl/2 _ 4| |„^|-7/2^2 12||5 |„^|-13/2^4 ^ q (^6)^ ^ 

(212) 

For X 00, Janson and Louchard [25| found asymptotic expansions begin- 
ning with 

Uo{x) - 3^3 xe-3-V2 (^i_l_^-2_A_^-^ + ,,,y (213) 
1 - f^n,e(a:) ~ V3 A _ ^ ^-2 ^ _L ^-4 ^ (214) 



23. Brovirnian double meander 

Define 

5dm(i) S(t) - min B(ii), te[0,l]. (215) 

0<ii<l 

This is a non-negative continuous stochastic process on [0, 1] that a.s. is at a 
unique point t e [0, 1] (the time of the minimum of B{t) on [0, 1]). It is well- 
known that T has an arcsine (= Beta(^, ^)) distribution with density 7r~i(t(l — 
t))~^/^; moreover, given r, the processes Bdm{T — t),t£ [0, r], and BdmiT + t), 
t G [0, 1— t] , are two independent Brownian meanders on the respective intervals, 



see Williams [57|, Denisov [13| and Bertoin, Pitman and Ruiz de Chavez [6|. 
(This may be seen as a limit of an elementary corresponding result for simple 
random walks.) We therefore call -Bdm a Brownian double meander. 
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Let Bd-m '■— Jq Bdinit) (it be the Brownian double meander area. By the 
definition, we have the formula 

Sdm:-^ Bd,n{t)dt = J^ B{t)dt- mm^B{t); (216) 



note the analogy with (fT6|) for Sox and Sbr- As in Section [3l there are further 
intere 
have 



interesting equivalent forms of this. Since B is symmetric, B = —B, we also 



Bdra = max B{t) - / B{t)At = max B°(t), (217) 

0<t<l Jq 0<t<l 

where ^ 

B°(t) := B{t) - I B{s) ds, < t < 1, (218) 
Jo 

is a continuous Gaussian process with mean and integral identically zero; its 
covariance function is, by straightforward calculation, given by 

Cov(B°(s),B"(t)) = i-max(s,t) + i(s2 + i2)^ s,ie[0,l]. (219) 

Let Vdm(<) := Ee"*^"^" be the Laplace transform of Bdm- Majumdar and 
Comtet [sll give the formula (using C(aj) for our r"^), where rj = AI(aj)/Ai'(aj) 
as in Section [22l 

C30 

i=i 

they then derive from it 

'e--Vdm(y2t3/2)dt= (^^)', a;>0, (221) 
using the partial fraction expansion, cf. (|2Q8[) . 

Mf)^ -V -V (222) 

(There are only quadratic terms in this partial fraction expansion, because 
AI'(z) = — Ai(2;) and Ai"(2:) = zAi(z) vanish at the zeros of Ai; hence there 
is no constant term in the expansion of AI(z)/Ai(z) at a pole aj.) Conversely, 
we will prove (|221l) by other methods in Section [271 and then p20p follows by 
P22p and a Laplace transform inversion. 

Majumdar and Comtet 35] found, by a Laplace transform inversion of (|220p . 
the density function fdmix) of Bdm as 

1/3 ^ 1 zL ^ A 

fdn.{x) = ^|a,|r|e-'^^ ([/(-,-;«,) +2[/(-^, (223) 

i=i 
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where U is the confluent hypergeometric function and Vj = 2|ajp/27a;^, cf. (|9ip 
and ([5^ . 

Asymptotic expansions of fdn\{x) and the distribution function i^dm(2;) as 
a; ^ follow as in [s^] from (|223[) . again using the asymptotic expansion of U 
[3^ . (9.12.3)]; only the term with j = 1 in (|223ll is significant and the first terms 
of the asymptotic expansions are, as a; ^ 0, 

(224) 
(225) 

For X — > cxD, Janson and Louchard [2^ found asymptotic expansions begin- 
ning with 

/d„.(x) ^ ^ (^1+1^-2^^,-4^;)^ (226) 



/3^ V 6 9 

The weaker statement 

- lnP(6dm >x)^- ln(l - Fdm{x)) ix^ /2, as a; ^ 00, (228) 

had earlier been proved by Majumdar and Comtet [35| using moment asymp- 
totics, see (I235P below. 

Define Wn by, in analogy with (|18ip . 



2-"/2n! 
r((3n + 2)/2) ' 



E-^dm = 7^77:;— ^^W/'n: "> 0. (229) 



An asymptotic expansion of the left hand side of (|22ip yields, arguing as for 

o, inni), (HUD and ([tosd . 



Al(a:) - ^ °" 
Ai'(a;) 



^(-l)'=VFfca;-3'=/2-i^ x^c5o. (230) 



By (IT55)) . (dnSl), PIM)1 and (fTM)) . this yields the equahty (for formal power 
series) 

i^s^m^) ' - f ' = (e q^A ' • (231) 



Consequently, 



W'n^^QjQn-j, n> 0, (232) 

J=0 
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dm \ ; dm 24 140 V TV "^^ 3840 



2 4 64 32 

Wq =1 = 12 W2 = 272 VK3 = 8856 W4 = 379264 

Table 5 

Some numerical values for the Brownian double meander. 



or, as given by Majumdar and Comtet [s^ (where IESJ,-,-j is denoted by /j,„), 
^^."n = ± (:>(^> fc|0±l) E^™ . (233) 

m=0 ^ ^ 

From (|232p and (|207p follows the asymptotic relation 

W^n ~ 2QoQn = 2g„ ~ 2V3(— j , (234) 
and thus by (|229p . see also Janson and Louchard [2E«] . 

EB2„~2V2(-) . (235) 

Some numerical values are given in Table O see also Majumdar and Comtet 
We define W* := 2^"Wn] these are integers by (i232| . 



24. Positive part of a Brownian bridge 

Let a;+ := max(x,0) and x_ ;= (— a;)+ and define 

/ Sbr(i)±di; (236) 
Jo 

thus is the average of the positive part and B^^ the average of the nega- 
tive part of a Brownian bridge. (Perman and Wellner 38| use Aq for B^^..) In 
particular, 

Sbr=6++Sbr. (237) 

and the difference is Gaussian: 



^br-^br= f' B^,{t)dt^N{0,l/U). 
Jo 



(238) 



By symmetry B^^. = B^^., so we concentrate on B^^.. Let V'br(^) ]Ee~*^br be its 
Laplace transform. Perman and Wellner [s^ gave (using the notation for 
•0^^.) the formula 

r e-^'i,l{V2t'/')^^2-—4^^^—r^, x> 0, (239) 
Jo a;i/2Ai(x) - Ai'(x) 
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see also Appendix IC.2I 
Define D+ by 



/¥2-"/2n! 



(Perman and Wellner [s^ use for our and /i+ for E(Sj^j.)".) An asymp- 
totic expansion of the left hand side of (|239p yields, arguing as for e.g. ([55]) and 

^^'^''^ - f;(-l)^-i?+:E-3'=/2-i/2, ^^oo, (241) 



x^/^k\{x) - Ai'(x) 



fc=0 



and thus by ^02]) . ^03]) . ^83]) . ^06]) . pOT]) . 

which leads to the recursion by Perman and Wellner [38| 

n _^ 

Using (frnS]) and (fn7|) . we obtain from ([M^ also 

E_j_ ^ 2 2 ^ ^ ^. Q Dj^z 

and the corresponding recursions 

71 

A— 1 

n 

= Dn~^Y. D^^n-k > 0- (246) 



fc=l 



Some numerical values are given in Table [6l see further Perman and Wellner 
111 (but beware of typos for fc = 5 and 7). We define D+* := 2^'^D+; these are 
213. 



integers by 
Tolmatz 



53l | gives the asymptotics, see also Janson and Louchard [25] and 



compare (|156p and (|155p . 

E{B+r ~ i EB^, ^ (— ) as n ^ oo, (247) 

or, equivalently, 

I':-iA.^i^(^)V-l)--^(|)" ^n^^. (248) 
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E(B+)0 = 1 EB+=i,/f E(B+)2 = i- E(B+)3 = ^/f e(B+)4 



2 20 4096 V 2 18480 

D+ = l D+ = l D+ = ^ D+ = ^ D+ = ^ 

" 1 8 2 32 3 5;^2 * 2048 

D+* = 1 D+* = 1 D+* = 6 0+* = 71 r>+* = 1266 

Table 6 

Some numerical values for the positive part of a Brownian bridge. 



Joint moments of B^^ and B^^^ can be computed by the same method from 
the double joint Laplace transform computed by Perman and Wellner (sl] , see 
((5551) in Appendix [C2] below. Taking A = 1 in (I^M)) . the left hand side can be 
expanded in an asymptotic double series, for ^, r/ \ 0, 

k,l=0 yJ^TTK.L. 

Using (|101[) and (|108[) , we find from (|369p the identity (for formal power series) , 
after replacing f by — ^ and 77 by —rj, 



k,l=0 



= 2 



Y.Zo<e . EZoC^'kV' V (250) 

k=0 "fe^ Ek=0 O^kV" 



1 



Note that r] = £, yields (fW)) ( using (|135p and (|152p ) and that rj — yields 

(EH. 

Define D^^ as the coefficient of ^'"'77' in the left hand side of (|250p . Thus 

E((B+f(B^V) = -4— ,\ Df,, kj>0 (251) 

vv br; V hru r((3(fc + z) + 1)/2) '^'^ ' - ^ ' 

and 

°° 1 1 

^fti^'^' = ^E^^WT^ ^ T^E^rWT^' ^^^^^ 

which yields the recursion formula, with Dqq = 1, for fc + I > 0, 



i=i j=i 
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We have = Z)*^ = D+ and 

n 

E^t.-^=^"' "^0- (254) 

fc=0 

We remark further that if n = 2m is even, then (|25ip and (|238p yield 

.X:(-ir-f")iEK)^M"- 



J=0 ^ j=0 



2"/2r((3n4 


-l)/2) 






2"/2r((3n-f- 


-l)/2) 


2"/2r((3n4 


-l)/2) 



¥n! 2"/2(„/2)! 



-Smo — m 



(6m)! 

(3m)! m!' (255) 

while the sum vanishes by symmetry if n is odd. Consequently, rj = — ^ in (j250p 
or ([252)1 yields 

V o-8mo-m_(M]_.2m ^ ^ ^ 

i-o (3m)! m!^ -2E™=o^2™e" 1 - 2 E™=i ^2™?2™ ' 

(256) 

which can be written as a recursion relation for i^2m with even indices only. 

Some numerical values are given in Table[71 In particular, as found by Perman 
and Wellner [38], 

Cov(Si.^..)-Tl)-li8' (2^^) 
Corr(e+ , 6,,) = « -0.636791. (258) 

We do not know of any formula for the density function /^^ ov distribution 
function Fj^, except the Laplace inversion formulas in Tolmatz [53] and Janson 
and Louchard ^25i] (which prove that /j^ exists and is continuous). For a; — > oo 
we have by Tolmatz [Hlj and Janson and Louchard [l^ asymptotic expansions 
beginning with 

\/6 _g3.2 A 1 _2 

-=e ^ IH X 

./^ V 36 5184 

l_F+(x)~^=x-ie-6^' (l-^x-^ + -^x-^ + ..\. (260) 
2\/6^ V 18 5184 J ^ ' 

We do not know any asymptotic results as x ^ 0, but the results on existence 
of negative moments in Section suggest that f^.^{x) x x~^/^. More precisely 
we conjecture, based on (|317p . that 

- (2^)3 (261) 
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br 
1,0 



1 



eb; 



3,0 



D, 



0,0 



D. 



3,0 



EBf' = — 



27r 

IF 

1 

2"0 

71v^ 
8192 



3 

32 
71 
512 



EB" 



EBi 



eb: 



EB: 



3,1 



V27r 
16 
1 

120 

13\/2^ 
24576 
1 

2880 

1 



Df, = — 



1 

3"2 

13 
512 

77 
2048 



EB" 



EBi 



EB^ 



EBf 



= — 



Table 7 



1 

20 

13v^ 
24576 

43 
332640 

17v^ 
1835008 
3_ 
32 
13 
512 
43 
2048 
255 
8192 



eb; 

EB^ 
EB^ 
EBj^ 
Di 
D 
D. 
D. 



br 
1,3 



± 
1,3 



± 
2,3 



± 
3,3 



71 V2^ 
8192 

1 

2880 

17v^ 
1835008 
11 

3564288 

71 
512 

77 
2048 

255 
8192 

3025 



65536 



Some numerical values for the positive and negative parts of a Brownian bridge. EB^^' is an 



abbreviation for ¥.{B'^^)*' {B^^f . 



25. Positive part of a Brownian motion 

Define ^ 

^bm / B{t)± di; (262) 
Jo 

thus is the average of the positive part and B^^^ the average of the negative 
part of a Brownian motion on [0, 1]. (Perman and Wellner [s^ use A+ = ^"""(l) 
for /Stluj-) In particular, 

Sbm = B+^ + B^^. (263) 
Note also that the difference is Gaussian: 

^b+m - ^bm = /' B[t) At ^ N{0, 1/3). (264) 
Jo 

Again, there is symmetry: B^^^ = B^^, and we concentrate on B^^. Let iphmi^) 
E e^*^bm be its Laplace transform. Perman and Wellner [s^ gave (using the no- 
tation 5*+ for i/'i^jjj the formula 

see also Appendix lC.il 

Define L+ by 



2-"/2n! 
r((3n + 2)/2)' 



= wT:;:-^^^™ , n> 0. (266) 
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= 1 


bm 


_ 72_ 




^ ''■"^ 1260 




6989 
40960 




= 1 




1 

~ 2 


2 16 


^ 64 


Li 


20967 
1024 




= 1 




= 4 


L+* = 68 


L+* = 2008 


Li* 


= 83868 



Table 8 

Some numerical values for the positive part of a Brownian motion. 



(Perman and Wellner [38] use for our and for El'Sbm)"-) asymp- 
totic expansion of the left hand side of (|265p yields, arguing as for e.g. ([55|) and 

The asymptotic expansions pM]) . p03)) . pl3)) of Ai(a;), Ai'(a;) and AI(x) yield, 
using (fTOe]) and (fTUT]) . the identity 

which leads to the recursion by Perman and Wellner [s^ 

n ^ 

L+ = ia„ + i/3„ + ^ ^^afcL+_„ n > 0. (269) 

k=l 

Some numerical values are given in Table [51 see further Perman and Wellner 
(ssj . We define L+* := 2^"'L+; these are integers as a consequence of (|298l) 
below. 

Janson and Louchard [i^ give the asymptotics, compare ()191|) and ()192p . 

1 /n\"/2 



E(6+J"^iEeC^^^(^-j , (270) 

4t ~ 5-^" ~ — j , as n ^ oo. (271) 

Joint moments of B^^^ and B^^^ can be computed by the same method from 
the double joint Laplace transform computed by Perman and Wellner jsl] , see 
(|36ip in Appendix IC.ll below. Taking A = 1 in (|36ip , the left hand side can be 
expanded in an asymptotic double series, for ^ \ 0, 

- E(-l)'^'^^K^^ + 0^((^^'-)'(^-")')^'^- (272) 
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Using (UnSl), HOSl), (dial), (UnSI), ^Bi, dUi, we and from ([361]) the 

identity (for formal power series), after replacing ^ by — ^ and rj by —77, 

2fc/2+i/2 .3(fc + ;) 



^00 



fc=o"fc^ J2k=o^krJ vEfc=o"fc^ L,k=o^krJ 

-2EZoKke~2j:z,Kkv'- 

Note that 77 = ^ yields pM]) and that 77 = yields dUS]). 

Define Lfi as the coefficient of ^'^jy' in the left hand side of (|273p . Thus 



T{3{k + l)/2+l 

and 



HiKji^a) - wt. , n/o fc'^ ^ (274) 



(275) 



kd=0 

which yields the recursion formula, with Lqq = 1, 
k I 



Lli = E ^^^U.i + E ^^^M-. + ¥i^Q^ + (276) 



where (5 is Kronecker's delta. 
We have L^g = L^^ = and 



T^^tn-k^Lr.^ «> 0. (277) 



fe=0 



We remark further that if n = 2m is even, then (|274p and (|264p yield 

^ 2"/^r(3n/2 + l) nl ^ (3m)! 

n! 2"/2(n/2)! 3"m!' (278) 

while the sum vanishes by symmetry if n is odd. Consequently, rj = — ^ in (|273p 
or ([275)1 yields 

(3m)! 2m ^ Sm^O Q2mC^" .„„„^ 
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EB°'° = 

bm 


1 


EB?'^ = 

bm 


V2 


EB°'^ = 

bm 


17 

96 


EB°'^ = 

bm 


251^2 
1260v^ 


EBi'" = 

bm 




ESi'^ = 

bm 


1 


EBi'^ = 

bm 


V2 


EB}'^ = 

bm 


149 




96 


3150F 


122880 


EB^'" = 

bm 


17 


ebI'^ = 

bm 




EB^'^ = 

bm 


109 


ebI'^ = 

bm 


31^2 


96 


3150F 


368640 


300300v^ 


ESi*'" = 

bm 


251v'2 




149 


EB?''^ = 

bm 


31v^ 


EB't'^ = 

bm 


9391 


126O0F 


122880 


300300v^ 


660602880 


r± _ 
-'"0,0 - 


1 


r± _ 
-'"0,1 - 


1 
2 


r± _ 
^0,2 - 


17 
16 


-'"0,3 - 


251 
IT 


r± _ 

^1,0 — 


1 
2 


lU = 


1 
8 


r± _ 
^1,2 — 


3 
16 


r± _ 

-'"1,3 — 


149 
256 


r± _ 


17 
16 


T± - 
^2,1 — 


3 
16 


r± - 
^2,2 — 


109 
512 


r± _ 
^2,3 - 


279 
512 


r± _ 

-^"3,0 - 


251 
IT 


T± - 
^3,1 - 


149 
256 


T± - 
^3,2 - 

Table 9 


279 
512 


r± _ 
^3,3 - 


9391 
8192 



Some numerical values for the positive and negative parts of a Brownian motion. EB^^^ is 
an abbreviation for E{B^^_^J'' ('^bm)' ■ 



a relation bet"ween Qj and Kk with even indices only. 

Some numerical values are given in Table[9] In particular, as found by Perman 
and Wellner 

Cov(S+,S,-J-^-^, (280) 
Corr(S+, - « -0.567185. (281) 

We do not know of any formula for the density function /j^^ or distribution 
function F^^, except the Laplace inversion formulas in Janson and Louchard 
(25j (which prove that exists and is continuous). For x ^ oo we have by 
Janson and Louchard [25] asymptotic expansions beginning with 



fUx) ^ e-'^'/' f 1 + _L ,-2 _ _^ ,-4 + . . . ] ^ (282) 



V 36 648 

We do not know any asymptotic results as a; ^ 0, but the results on existence 
of negative moments in Section suggest that f^^{x) x x~^^^. More precisely 
we conjecture, based on (|319|) . that 

f+Jx) I 2^5/^3l/6^-3/2^(l/3)a;-2/^ x ^ 0. (284) 
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26. Convolutions and generating functions 

Many of the formulas above can be written as convolutions of sequences, as is 
done by Perman and Wellner [s^. Tiiis gives tiie simple formulas below, where 
a letter X stands for the sequence (X„)[f and X *Y is the sequence defined 
by {X * Y)n ■= J2k=o-^kyn-k- Alternatively, the formulas can be interpreted 
as identities for formal power series, if we instead interpret X as the generating 
function J2'^=o -^nz" and * as ordinary multiplication. We let 1 denote the 
sequence (Sno) — (1, 0, 0, ... ) with the generating function 1. 

Note that the sequences denoted by various letters below, except K, all have 
Xq = 1, and thus the generating functions have constant term 1; the exception 
is Kq = —\: which explains why —2K occurs frequently. (It would be more 
consistent to give a new name to ~2Kn and use it instead of Kn in our formulas 
(mutatis mutandis), but we have kept Kn for historical reasons, and because all 
numbers Kn with n > 1 are positive. See also pill) , which shows that {6n—2)Kn 
is more natural for some other purposes.) 

In the following list of formulas, we also give references to the corresponding 
equations above. Many of these formulas can also be found in Perman and 
Wellner 





-2a *K 


= a' 


(|1()8|), dZOl) 


(285) 




a' *D 


= a 


(|137l). (11481) 


(286) 




-2K^D 


= 1 


(|152l).p53p. (11541) 


(287) 




a' * L 


= 


(I185D. (I186D 


(288) 




a* Q 




(I196D. (I201D 


(289) 




-2K^L 




(|202l). (12031) 


(290) 




D*Q 


= L 


(12041). (12051) 


(291) 




Q*Q 


= W 


(I231D. (12321) 


(292) 




a * a *W 


= /3*/3 


(|23ip 


(293) 


(a 


■fa') 


= 2a 


(I242D. (I243D 


(294) 


D+- 


-2K*D+ 


= 21 


(12441). (12451) 


(295) 


D+ 


+ D*D+ 


= 2D 


(I244D. (I246D 


(296) 


{a 


+ a')^L+ 


= a + 13 


(I268D. (I269D 


(297) 


D+ 


*Q + D+ 


= 2L+ 




(298) 



The final relation ((298)) follows from ((294|) . (|289l) and ((297)) by simple algebra 
with generating functions. 

27. Randomly stopped Brownian motion 



Several of the formulas above have simple interpretations in terms of a Brownian 
motion stopped at a random time, which has been found and used by Perman 
and Wellner [38 1. 
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Let Sbm(i) '■— /q \B{s)\ ds; by Brownian scaling, Bhra{t) = i^^^^bm- Let fur- 
ther Z ^ Exp(l) be an exponentially distributed random variable independent 
of the Brownian motion B. The area of the (reflected) Brownian motion stopped 
at Z is thus 

Sbm(^) = ^'/'Sbn. (299) 

with moments given by, see (|18ip . 

E(;Bbm(^))" = EZ3"/2e6[,'^ = 2-"/2„!L„. (300) 

Thus, 

E(V2Bbm(^))" , , 

Ln = ^ , , n>0, 301 

and the generating function X]fc!Lo ^fc-^*^ (|185p is the exponential generating 
function of the moments of \/2Bh■al{Z)^ i.e. the moment generating function of 
\/2 -Bbm interpreted as a formal power series. (This series diverges, so the moment 
generating function in the ordinary sense Eexp(z\/2Sbm(-^)) is infinite for all 
z > 0.) Furthermore, the Laplace transform of \/2 BhraiZ) is given by, see p62l) 
below, for ^ > 0, 

^,-,V-2....iz, ^ J- e-Vb4V2et3/-) d^ = ,(P/3) - (302) 

which is equivalent to (|173p by a change of variables. 

Similarly, with obvious definitions, (j266p can be written 

, E(V26+ (Z))" , , 

and the generating function X^fe^o ™ (|268p is the moment generating func- 

tion of v^'Bbni interpreted as a formal power series. The Laplace transform of 
V2B^^{Z) is given by, see ([363]) below, for ^ > 0, 

Ee-«^«.^»(-)- re-v+ (V2a^/^)d^-M::!^!l±i:^^!M:!^ 

"io '^'^-('^^^^ ^ " Ai(C-2/3) - Cl/3Ai'(e-2/3) ' 

(304) 



which is equivalent to ()265|) by a change of variables. 

Furthermore, if Bhr{t), ^^A^) ^^'^ Bhm{t) are defined in the obvious way 
using a Brownian bridge and a Brownian meander on the interval [0,t], and 
Z' - r(l/2), then 11351), ^Mi) and can be written (since E(Z')^"/^ = 

r(3n/2 + l/2)/r(l/2)) 

E(V26br(^'))" . ^ 

Dn^^ , " , n>0, 305 

, E(\/2S+,(^'))" , ^ 

At = , ' , n > 0, (306) 
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Qn = ^ i n > 0, (307) 

and thus Y^k^kz^, Y.kDtz'", J2kQkz'" in ^37^, IMl, (Uli) are the mo- 
ment generating functions of y/2Bhr{Z'), V^B^j.{Z') and ^/2Bmc{Z'), respec- 
tively, interpreted as (divergent) formal power series. Since Z' has the density 
^-i/2^-i/2g-t^ i > 0, the Laplace transforms of these variables are given by, 
using (P7D1) . (P7T|) . (|5551) below, for ,^ > 0, 

Ee-«^«-(^') = e-V'b.(V2Ci3/^)-|= = -^'/'^^^y (308) 

di Ai(r^/^) 



Ai(e-2/3)_^l/3Ai'(^-2/3)' 

(309) 



IEe-«^«'-(^')=^°°e-Vmc(V2et^/^) J= = r^/^§|p^, (310) 

which are equivalent to (fT33)) and 1134)) . (|239)) and p93)) . 

The reason for the introduction of Z' is that, as is well-known and explained 
in greater detail by Perman and Wellner [s^, if gz ■— max{i < Z : B{t) = 
0}, the last zero of B before Z , then conditioned on gz, the restrictions of B 
to the intervals [0,32] and [gz,Z] form a Brownian bridge and a Brownian 
meander on these intervals, and the lengths gz and Z — gz of these intervals 
both have the distribution of Z' . Moreover, the restrictions to these two intervals 
are independent (also unconditionally, with gz and Z — gz independent), and 
thus BhmiZ) = Bhiigz) + Bme{Z — gz), where the two variables on the right 
hand side are independent with the distributions of Bhr{Z') and BjneiZ'), which 
immediately yiel ds (|29ip . Considering only the positive part of B, we similarly 
obtain ((298)) 

For Bex, the slightly different (O can be written 

E(V2Bcx{Z'))" , ^ 

(on — 2) n\ 

Hence, the proper analogue of L„, i+, £>„, -D+ and Qn is (6n — 2)7^„. Note 
that it follows from (|10ip that, as x — > 00, 

y (_l)n(6n - 2)if„X-V2-3n/2 ^ ( ^ ^ ( ' _ ^ 

(312) 

the Laplace transform appearing in ([M|) : further, in analogy with p08p - (|310p . 
([M]) is equivalent to, for ^ > 0, 

Ee-«^^-(^') = e-Vex(V2Ci^/2) J= = 2r^/^ ( 1^(^2/3) )' " 

(313) 
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For Bdin we similarly let B^^^ be a Brownian double meander on [0,t\. Then 
^dm ^-^^ ^ unique zero at a time /i = tZ where r ^ Bcta(i, i) is indepen- 
dent of Z; hence fj, and Z — \i = (1 — t)Z are two independent random variables 
with a r(l/2) distribution. Moreover, see Section [23l given /x and Z — \i, 
consists of two independent Brownian meanders on the intervals [0, /u] and [/x, Z], 
with the first one reversed. Consequently, the Brownian double meander area 
Bdm{Z) := (^) equals the sum of two independent Brownian meander 

areas, both having the same distribution as BmdZ'). Hence, by (|310p . 



-2/3^ 



VAi(e-2/3) 

which is equivalent to (|221|) by a change of variables. Furthermore, 



(314) 



W„ = , " , n > 0, 315 

n! 

and the generating function X^fe^o^'^-^'^ (|23ip is the moment generating 
function of V^Bdm interpreted as a (divergent) formal power series. 



28. A comparison 

For easy reference, we collect in Table [TU] the first two moments of the various 
Brownian areas treated above, together with the scale invariant ratio of the 
second moment and the square of the first. 

Note that the variables B^^, Bdmi Bj^a Bi^m have quite small variances; 
if these variables are normalized to have means 1, their variances range from 
0.061 to 0.325 (in this order, for which we see no intuitive reason), which means 
that these variables, and in particular Bex, typically do not vary much from 
their means. 



29. Negative moments 

We gave in Section [18] formulas for negative moments of Bex due to Flajolet 
and Louchard 17]. These can be generalized using the results in Appendix [Bl 
First, ([339]) and (IMO]) yield by (gl, (fTSS]) . (fTTS]) . pM]) . (|22T|) . (1239]) and (|265l) 
expressions for the negative moments of order an odd multiple of — ^ for Sex, 
Bhi, Bmc, B^^, and a multiple of — | for Bhm, Bdm, ^hin- ^'^^ example, this yields 
the values in Table [TT] (including a repetition of (|124p and p25p ): note that 
p26p enables us to rewrite the values in Table [TT] in other forms. 

More generally, by ([34T]l and ([84]), ([133]) . ([T73]) . ([T93| . ([22T]) . where the right 
hand side is analytic at 0, the random variables Bex, Bhr, Bhm, ^mo, Bdm have 
negative moments of all orders; thus, ^B^^, ^B^^., MB^^, EB^^, ^'^dm 
entire functions of s. 
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X 


EX 




EXV(IEX)2 


Bex 


4 


5 

12 


10 

— Ki 1.06103 

37r 




8 


7 

60 


56 

Ri 1.18836 

IStt 


^bm 


2v^ 
30F 


3 
8 


277r 

Si 1.32536 

64 


Bmc 


3v^ 

8 


59 
60 


472 

Si 1.11291 

1357r 


Sdm 




17 
24 


177r 

Si 1.11265 

48 




3v^ 


17 
96 


5l7r 

Si 2.50346 

64 


br 


16 


1 

20 


32 

— Si 2.03718 
57r 



Table 10 

The first two moments for seven Brownian areas. 



For Sj^j. we see that the right hand side of (|239p is finite at x = 0, but that 
its derivative is ~ cx"^/^ as a; ^ for some constant c 7^ 0. Hence, (I34ip 
with m — 1^ V — 1/2 and '^{x) equal to the right hand side of (|239p . shows 
that for < s < 3/2, the moment E(6+ is finite if and only s > 1/2. 
Consequently, the negative moment E(Sj^j.)~^ is finite if and only if s < 2/3, 
and E(;8^j,)'' is an analytic function in the half-plane Re s > — |. More precisely, 
a differentiation of (|239p shows that, for this 5", 

-*'(x)=0F(^)Vv2 + O(i), :,^o, (316) 

and thus the left hand side of (|34ip with m — 1 has the residue -y/7r(Ai(0)/Ai'(0))^ 
at s = 1/2, and consequently (|34ip yields 

p.. vtn+Y fAm\' ^ 2V33-^/^r(l/3)^ _ 2V335/6r(l/3)5 

^=-3 ^ i^'-^ ~ r(2/3) Ui'(0)J r(2/3)3 (27r)3 

(317) 

For we see that the right hand side of (|265p is ~ cx as a; — > for 
some constant c > 0. Hence, p38p with v — 1 and ^'(x) equal to the right hand 
side of shows that for < s < 1, the moment E(S+^)-2(i-'*)/3 ig finite if 

and only s > 1/2. Consequently, the negative moment ¥,{B^^-/)^^ is finite if and 
only if s < 1/3, and ^{B^^^ is an analytic function in the half-plane Res > — |. 
(This also follows from the considerations in Section [^H which imply that, for 
s > 0, E(6+^(Z))-'^ is finite if and only if E(S+ (Z'))"" = HKX' HZ')'"^"'^ 
is, and E(Z')"3s/2 < jf ^nd only if s < 1/3.) More precisely, (|265p shows 
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21/63^ /Ai'(0)\2 

r(i/3) V Ai(o) ; 


= 2" 


-17/63l9/6^-5/2p(2/3)5 , 


» 1.184 


br 


2-5/63^ Ai(0) 
r(l/3) -Ai'(O) 


= 2" 


-ll/637/6^-l/2r(i/3) , 


s 1.528 


]gg-l/3 
^ '-'me 


2-«/«3Vi AI(0) 
r(l/3) Ai(0) 


= 2" 


-ll/637/6^-l/2p(2/3)2 , 


» 1.046 


E(S+)-V3 


2i/«3V¥ Ai(0) 
r(l/3) -Ai'(O) 


= 2" 


-5/637/6^-1/2^(^/3) ^ 


a 3.056 


Esr^/^ 

bm 


2-2/33 AI(0) 
r(2/3) -Ai'(O) 


= 2" 


-2/33I/3 '^(1/'^) 

r(2/3) 


» 1.797 


EST^/^ 

dm 


2-2/33 /AI(0)y 
r(2/3) lAi(0)j 


= 2" 


-2/33l/3p(2/3) r 


s 1.230 




= 00 









„-i 35/2r(2/3)6\ 

ESex = 3V^ ( 1 J^3^ ) 

ESj^/ =2-1/237^ 

E.- ^2-V23Vi(l-^^^!§HZ^) 

E(i3+ )-i = 00 



Ee-*/3 

bm 


= 2-4/33i3/6^-ir(l/3) 


» 3.658 


ES,-*/^ 

dm 


- 2-i/33i3/67r-ir(2/3)2 (1 3i/2r(2/3)3 ^ 


« 1.580 


Esr2 

bm 


3 , 3i/2r(l/3)3 
2n 


« 8.300 


EB72 

dm 


g 35/2r(2/3)3 27r(2/3)6 
TT 27r2 


« 2.112 


Table 11 



Some negative fractional moments. 



« 1.693 
« 3.760 
« 1.186 
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that, for this 

*(^)=-^.^"'/'+0(l), x^G, (318) 

and thus the left hand side of p38|) has the residue — Ai(0)/Ai'(0) at s = 1/2, 
and consequently (|338p yields 

Res 1 mB+ Y = ^ ' = 

^es,__ yo^^, r(l/2)r(l/3) -Ai'(O) V^r(2/3) (319) 
= 2-5/63i/6^-3/2r(i/3). 

Furthermore, the right hand sides of ()239|) and (I265P have expansions in 
powers of x^/"^ at a; = (beginning with x~^/'^ for (|265p ). Hence, in p4ip with 
\l/(a;) equal to one of these functions, the left hand side extends to a meromorphic 
function in < Re s < m+u with possible poles (all simple) only at ^, |, . . . , and 
thus the same holds for the negative moment EX~^^'"+'^~*^/'^ on the right hand 
side of (|34ip . Since m is an arbitrary non-negative integer and v = 1/2 for B^^ 
while I' = 1 for B^^, it follows that '&{B^^Y extends to a meromorphic function 
in the complex plane with poles only at — |, — |, —2, i.e., multiples of 
— |, and that ¥,{B'^^^'' extends to a meromorphic function in the complex plane 
with poles only at — i, —1, — |, i.e., odd multiples of — i. Furthermore, 
all poles are simple. (We cannot prove that all these points really are poles; 
it is conceivable that some actually are regular points because some terms in 
the expansions of "^{x) might vanish due to unexpexted cancellations, but we 
believe that all actually are poles. At least, — | and — ^ are poles by the results 
above.) 

Recall also the relation (|122p in Section [18] between moments of Bex and 
the root zeta function A of Ai. For the Brownian bridge, there is an analogous 
relation between moments of Sbr and the root zeta function of Ai' defined by, 
cf. (fT20l) . 

OO 

A(s) ^ |a^.|"^ Res > 3/2. (320) 

Since \a'^\ - \aj\ - {i-Kj/if/^ P, 10.4.95] (in fact, \aj_i\ < \a'^\ < \aj\ for j > 2), 
the sum converges and A(s) is analytic for Res > 3/2. By (|335p in Appendix IB] 
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and (|159p . for Res > and with absolutely convergent sums and integrals, 

/>oo 



= 2-1/6 / V|a;riexp(-2-i/3|a;|t2/3)^ 

= 2-^/«3^/^^|a;rW u(3^+i)/2exp(-2-i/3|a;|H)— (321) 

OO Q I 1 

= 2^/2-13^^ |a^.|-(3.+3)/2p|^£l_Li^ 

Since KB^^ is an entire function, this shows that A extends to a meromorphic 
function in the complex plane with, similarly to (|122p , 



Since the left hand side has no poles and no real zeros, it follows, arguing as in 

3 _3 _9 

2 ' 2 ' 2 ' 



Section [m that A(z) has simple poles atz=|, — |,— |, and zeros at 1, 



— 1, —2, —4, —5, —7, . . . , and nowhere else. 
It follows from (fT551) that, for |z| < |ai|. 



Ai(^) 
Ai'(z) 



^(-l)™+iA(m + 2)0™, (323) 



and thus the values of A(s) at positive integers s = 2,3,... can be com- 
puted from the Taylor series of Ai. (In particular, A(2) = — Ai(0)/Ai'(0) = 
3-i/3r(l/3)/r(2/3) and A(3) = 1. Furthermore, as just said, A(l) = 0.) This 
and (|322p give again the same explicit formulas that we have obtained from 
(|340p for the negative moments B^^^ when s is an odd multiple of 1/3. 

Note also that by (|133p . p38p and p2ip . we have the Mellin transform, for 
< Res < 1/2, 

-r(s)r(l-s)A(2-s) (324) 
= -- ^A(2-s), 



which also follows from ([T58l) : cf. (fT27l) and (fT28| . 

For Brownian motion, Brownian meander and Brownian double meander, 
more complicated analogues of the root zeta function A(s) appear in the same 
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way: for Brownian niotfon, using (|170p . X^jli ^jl^jl ^o^' Brownian mean- 
der, using (|209[) . ''ikil f^'' Brownian double meander, using (|220[) . 
X^jli l '^^ We do not pursue this further. 



Appendix A: The integrated Airy function 

We define, as in ([78]), using Ai(x) da; = f/3 0, 10.4.82], 



+ 00 



1 



Al(;z) := / ki{t) dt = - - / ki{t) di; (325) 

this is well-defined for all complex z and yields an entire function provided the 
first integral is taken along, for example, a path that eventually follows the 
positive real axis to -foo. Note that Al'(z) = —Ai{z) and that A1(0) — 1/3. 



Along the real axis we have the limits, by [ij, 10.4.82-83], 

lim AI(.t) = 0, (326) 

X — >-\-OQ 

/oo 
Ai(a;)dx = l. (327) 
-oo 

In terms of the functions Gi and Hi defined in [ij, we have, see pj, 10.4.47-48], 

AI(z) = 7r(Ai(z)Gi'(z) - Ai'(z)Gi(z)) (328) 
= 1 + 7r(Ai'(z)Hi(z) - Ai(z)Hi'(z)) . (329) 

Repeated integrations by parts give 

p+co r+oo 1 

Al{z)^ Ai{w)dw= -Ai"{w)dw 



w 



Ai'(z) f+°° 1 
Ai'(z) Ai{z) f+°° 2 



z 

f + C30 



r+oc 2 

/ — ^Ai(u;) duu 
Jz w-^ 



(330) 



Ai'(z) Ai(z) /•+°° 2 ^.„, ,^ 
^ ' w , ; Ai (u;)dw 



if we assume ] arg(z)] < n — 6 for some 6 > 0, we can here integrate for example 
along the horizontal line {z + t : t > 0}, and it follows from the asymptotic 
expansions (|102p and (|103p . which are valid for ] arg(z)] < tt — S [1, 10.4.59 and 
10.4.61] that there is a similar asymptotic expansion, as \z\ oo, 

oo 

AI(z) ^ -J-z-3/V2^''^/3 £(-l)^-/3,z-3'=/2, (331) 



k=0 



valid for ] arg(z) \ < tt — S and thus extending (|183p (which is given in Takacs 



y(z 

l49|; see also [l|, 10.4.82]). The coefficients f3k can be found by this procedure. 
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but it is easier to observe that we can formally differentiate (|33ip and obtain 
an asymptotic expansion of — Ai(z), and a comparison with (jl02p yields the 
recursion relation (|184p . See Table [3] for a few numerical values. 

Similarly, along the negative real axis, and more generally for —z where 
I arg(z)| < 27r/3 — 5 for some ^ > 0, there are asymptotic expansions [1, 10.4.60 
and 10.4.62] 

Ai(-.) ^ .-V. ( sin(^z^/^ + J) Y.{~lfc.,,z-''-"' 

- cos(^z'/' + J) Y.{-ira,,^,z-''^-yA^ (332) 

Ai'(-z) ^ -vr-V2 ( cos(|z3/^ + J) E(-l)'=4,z-^'^+i/^ 

^ fc=0 

9 °° \ 

+ Sin(2 3/2 ^ ^(_^).^,^^^^-3.-5/4y (333) 

fc=0 ^ 

which using (I330|) lead to, as \z\ ^ oo with | arg(z)| < 2tt/3 — S, 

/ 9 °^ 
AI(-z) ^ 1 - 7r-V2 ( cos(|z3/2 + J) 5:(-l)'=/?2feZ^^'=-^/^ 

^ fc=0 

9 °° \ 

+ sin(|z3/2 + J) ^(-l)'=^2.+iz-3'=-^/M ; (334) 

k=0 ^ 

given (at least along the negative real axis) by Takacs 'i^, (17)] (using his 
hk = (2/3)'=/3fe); see also 0, 10.4.83]. 

Appendix B: The Mellin transform of a Laplace transform 

For any positive random variable X with Laplace transform ip{t) := Ee^*^, 
and any s > 0, Fubini's theorem yields 

/>oo />oo 

/ t''-^i;{t)dt:^E f-^e-*^ dt^T{s)EX-''. (335) 

^0 "'0 

If one side is finite for s — sq > 0, and thus both sides are, then EX^" < oo 
for s — So and thus for < s < so; hence both sides of (|335p are finite and 
(|335p holds for all s with < s < sq, and more generally for all complex s with 
< Re(s) < So, with the terms in p35p analytic in the strip < Re(s) < sq. 

Lemma B.l. The following are equivalent: 

(i) X has negative moments of all orders, i.e., E AT^'* < oo for all s > 0; 
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(ii) t'-'^i:{t) dt<oo for all s > 0; 

(iii) V(0 = 0{t-^) as t ^ oo for all N. 

If these hold and further X has moments of all (positive) orders, then EX~* 
is an entire function of s; further, (|335p holds for Res > and the right hand 
side defines a meromorphic extension of t''~^'ijj{t) dt to the complex plane. 

Proof (i) ^ (ii) by (15551) . 

Since V(0 < V'(O) = 1, (iii) =^ (ii). 

Conversely, since -0 is decreasing, (ii) implies for every m > 

/■oo /'U 

oo> t'^^i^it) dt > ij{u) t''-^dt = s-\'i;{u), (336) 
Jo Jo 

and thus ip{u) = 0{u'~^). 

The last statement is obvious. □ 

Several formulas above, viz. (fTSS]) . (fT73ll . (fTM]) . ([22T|) . (|239| . ([MS]) . give 
a formula for a double Laplace transform of the form 

e-^*V(V2t3/2)t'^-idi = *(x), x> 0, (337) 

where = 1 or 1/2; see also Section [27] for an explanation of this form. (We 
have 1/ = 1 for Bhm and and = 1/2 for the others.) We then can argue as 
above twice (using ()335p the second time) and obtain for < s < i^, and if the 
result then is finite also for complex s with < Re s < i>, 

/>oo 

^■^{x)dx^T{s) / t-'i}{V2t^^^)t''-^dt 



^T{s)l I 2-(''-^)/3^2(.-.)/3-i^^(^)j^ (338) 



3 7o 

= ^^:^r(.)r(fcil)EX-(-)/3. 

Furthermore, a; = in ()337|) similarly yields, with ^'(0) := limj;^o ^(a;) < oo. 



*(0) ^ / V'(V2i=^/2^i''-Mi 

= 1 / 2"'/\'''/^''i;{u)du (339) 
3 Jo 

= r(— ) EX-^'^/s 

3 V 3 / 



which can be regarded as a limiting case of (j338[) for s ^ 0. In particular, 
4'(0) < oo if and only if E.X^^'^^^ < oo (which holds for all Brownian areas 
treated here except B^^). 
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More generally, by first taking the mith derivative of (|337p . we obtain for 
every integer m > 0, by letting a; — > and with 5'^™)(0) := linix^o 5'^™-' (a;), 



. w (340) 

and, for < s < m + z/ (and for < Re 5 < m + if the result is finite), 
(-1)™/ a;"-i^'(™)(x)da; 

Appendix C: Feynman-Kac formulas 

Many of the results above are derived from formulas obtained by simple ap- 
plications of the Feynman-Kac formula, following Kac 13; Slj . This is by now 
standard, but for completeness we give here a unified treatment of all the Brow- 
nian functionals discussed above. We basically follow the proof by Shepp 42] for 
the Brownian bridge. See also Jeanblanc, Pitman and Yor 27[ for a treatment 
of several related Brownian functionals. 

One version of the Feynman-Kac formula is the following. Let Bx{t) denote 
the Brownian motion started at x; thus Bx{t) = x + B{t)] we continue to use 
B{t) for Boit). 

Lemma C.l. Let V{x) > be a non-negative continuous function on (— cx),oo), 
let A > 0, and let (p^ and (j)- be solutions of the differential equation 

(l)"{x)^{V{x)+\)(j){x) (342) 



1a," 

2 



such that is bounded on [0,oo) and is bounded on (— oo,0]. Let w := 
(/)+(0)(/)'_(0) — (/)_(0)(/)^(0) and assume that w ^ Q. Moreover, let f be any 
bounded function on M. Then, for all real x. 



j e-^*E(e"io ^(^^(^»'^V(S,(i))) At 



= ^ (^M^) J <l^-{y)f{v) + 0-(a:) ^ 0+(y)/(y) dyj . (343) 

Proof. Although this is well-known, we give a proof for completeness. 
Note that if we define the Wronskian 

W{x) := c^+{x)c^'_(x) - c^^{x)c^'+(x), (344) 

then p42p implies that W'{x) = and thus W{x) is constant; hence 

W{x) = w. (345) 
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Assume first that / is continuous and has compact support, and define 

/X POO 
</)_(y)/(y)dy + 0_(a:) / dy. (346) 

Then G C^(M), and differentiation yields 

cj)^{y)f{y)dy + ^'_{x) / 0+(2/)/(y) dy. (347) 

OO J X 

Hence also (/)' G C^, and thus G C^. A second differentiation yields, using 

(j)"{x) = 2(F(x) + A)(/)(a;) - VK(a;)/(a;) 2{V{x) + A)(/)(a;) - u;/(a;). (348) 

Moreover, since / has compact support, ()346|) shows that if x is large enough, 
then 0(x) = (/ 0_/)(/)+(x), and thus (/)(a;) is bounded for large x; similarly (f){x) 
is bounded for x — > —oo. Consequently, </> is bounded on M. 

Let ^(i) and Y{t) be the stochastic processes X{t) At + /oT/(B,(s)) ds 
and ^ 

F(t) e-^(*)0(i3,(t)) + ^ / e-^(")/(i?,(7.)) dw. (349) 

^ Jo 

A straightforward application of Ito's formula, using (|348p . shows that 

dY{t) = e-^(*V'(S.(0) dB,(t), (350) 

and thus Y{t) is a local martingale. Moreover, since Xt > Xt and and / are 
bounded, Y{t) is uniformly bounded, and thus a bounded martingale. Hence, 

(t){x) =Y{0) ^EY{oo), (351) 

where, recalling that (f) is bounded and X{t) > Xt, 

/•oo 

r(oo) := lim y(t) - ^ / e-^(")/(S.(^))dw. (352) 

The result p43p now follows by (|35ip . p52p . and p46p . under our assumption 
that / is continuous with compact support. 

By a monotone class theorem, for example [2^ Theorem A.l], p43p remains 
true for all bounded measurable / with support in a finite interval [— ^, A], for 
any A < oo. Fixing an x with 0_(a::) 7^ and letting f{y) :— sign(0-|_(y))l[a:: < 
y < A], we see by letting A ^ 00 that \(j)+{y)\dy < 00. Hence is 
integrable over each interval (x,cxd), and similarly (jj- is integrable over each 
interval (— oo,x). It follows by dominated convergence, considering /(x)l[|a;| < 
n] — > f{x), that (j343p holds for all bounded measurable / on M. □ 

Remark C.2. Majumdar and Comtet [35| have given "physical proofs" of (|80p 
and (|209p using path integral techniques. This method is closely related to 
the Feynman-Kac method as formulated in Lemma FCli which can be seen as 
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follows, where we argue formally and ignore giving proper technical conditions 
on V and / for the validity of the argument. 

Let H be the differential operator Hcf) :— —^(t/' + Vcj). (This is an unbounded 
positive self-adjoint operator in L^(R).) Then (|348p can be written 

w 

H4, + \(jy=-f. (353) 

Hence, the right hand side of (|343p . which is :5<^(^) P46p . equals {H + 
X)~^f{x). Further, {H + X)~^f = e~'^*e~*^/ dt, so by the uniqueness theo- 
rem for Laplace transforms, (I343|) is equivalent to 

E(e-/o ^(^^(^»'^7(i?,(t))) = e-'"f{x), (354) 

which essentially is the path integral formula used by Majumdar and Comtet 
[35I I in their proof. 

We let ^, ry > and choose 

n^):=|^^r' (355) 
^ ' \V2v\x\, x<0. ^ ' 

It is easily checked that Ai(^/2^^/^a; -I- ^"^/^A) solves (p42l) on (0, cxo) and is 
bounded there; thus this can be extended to a solution (/)_|_ on R that is bounded 
on [0,00). Similarly, Ai(^—^/2 rj^^^x + rj^^^^X^ on (— oo,0) can be extended to 
a solution that is bounded on (— oo,0]. We thus have (p± as in Lemma FC. II 
with 

(j)+{x) = Ai{V2^^/^x + C^/^X), x> 0, (356) 

(j)^{x) = Ai{V2T]^/^\x\+'n'^/^X), x< 0. (357) 

As a consequence, 

w = -y2ryi/3Ai'(,^-2/3A)Ai(r'/'A) - ^e^''Ai\r^/'X)Ai{7]-^/'X). (358) 
We will now apply Lemma FC. II with this V to the different cases. 

C.l. Brownian motion 

Recall the definition (I262p of B^^^ and define the joint Laplace transform, for 

V'bm(e-'?):=IEe-««-,-')6,-„. (359) 

Thus Vb„.(0 = ^bm(^.0 and ^^(0 = ^^(^,0). 

Assume that ^, 77 > and use Lemma FC. II with V as in ([355]), 4>± as in ([356]) 
and ()357|) . / = 1 and a; = 0. By (|355p and homogeneity, 

/* V{B{S)) ds 4 t3/2 r V{B{S)) ds - V2t3/2eS+ , + 72*3/2^^-^^ (ggg) 

Jo -'0 
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and thus piS]) yields, recalling (|359l) and using (|356l) . ([BST]) . ([358]) and ([78l) . 

^ (^0+ (0) y ° ^- {y)dy + 0_ (0) y 0+ (y) dy^ 



_ / AI(r'/'A) AI(77-^/3A) A / ^V3Ai'(g-2/3^) ,^l/3Ai'(r^"2/3;^) X -1 

~ UV3Ai(5-2/3A) + 77i/3Ai(r7-2/3A) j Ai(e-2/3A) Ai(?7-2/3A) J 

(361) 

In particular, the special case 77 = ^ yields 

Similarly, letting 77 ^ we obtain, since AI(a;)/Ai(a;) ~ x^^/^ and Ai'(x)/Ai(a;) ~ 
as a: CX3 by (fTOS]) . (fTOS]) and (fT83l) . 

^ / AI(r^/^A) i/,A / gV3Ai^(^-2/3^) X -1 

Ui/3Ai(e-2/3A) ^ ; 1, Ai(r2/3A) + ; 

_ 1/3 AI(r'/'A) + A-V^gV3Ai(r'/'A) ^ ^ 

" ^ -ei/3Ai'(^-2/3A) + Ai/2Ai(^-2/3A) 

_ 1 Ai(r'/^A) + r'/^Ai/^AI(r'/^A) 
Ai(e-2/3A) ^ e'/3A-i/2Ai'(e-2/3A) ■ 

Taking,^ = 1 yields ((2651) . 
C.2. Brownian bridge 

Let Bx}y denote the Brownian bridge on [0, i] with boundary values Bx\{Q) = 
X and Bx}y{t) — y; it can be defined by conditioning B^ on {Bx{t) — ?/}, 
interpreted in the usual way (with a distribution that is a continuous function 
of y) . The standard Brownian bridge _Bbr equals i?g^Q ■ 

The left hand side of p43p can be written, by conditioning on Bx{t) inside 
the integral, 

e-At r Ee-/o'^(^a(^))d^/(y),-(i/-)V(2t)^di (364) 
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while the right hand side can be written, with xM y := max(a;, y) and x A y := 
mm{x,y), 

2 

(t)+{xyy)(j).{xAv)f{v)dv. (365) 



w 

Hence, using also Fubini in (I364p . 

Te-^'Ee^/o ^(^-(^"'^^e-f^-^)^/^^*)^ = -^xyy^^ixAy) (366) 
Jo \/2iTt w 

for a.e. y. However, both sides of (|366p are continuous functions of y (the left 
hand side by dominated convergence because Bx y{u) = B^^\{u) + yu j t in C[0, t] 

and thus Bx y^^ — > Bi% in C[0, t] as j/„ — > y). Hence, p66l) holds for every 
y€R. 

Take x — y — and F as in (|355p . Then, by homogeneity, recalling (|236|) . 
/* y(S^:i(s)) d. 4 t3/2 r as = V2t'/'^Bi^. + V2t'^'vB^,. (367) 

JQ Jo 

Define the joint Laplace transform, 

thus ^br(0 = ^iiLO and i^+iO = V'br(C,0). Then, §661) and jMD yield, 
using (i356l) - ([358ll . for C,?? > 0, 



(369) 



'e-^V^,(V2tn, V2t3/^^)-^ = -0+(O)0-(O) 

_ ^ / e'/'Ai'(r'/'A) 77V3Ai'(r;-V3A) \'^ 
Ai(C-2/3A) Ai(r/-2/3A) j ■ 

In particular, taking ^ = 77 > 0, 



which by simple changes of variables is equivalent to ()133|) (A = a;, ^ = 1) and 
(fT34l) (A = X, C = w/V2). Further more, letting 77 wc obtain 

dt _,,_i/2 Ai(r^/3A) 



/ e-^*^+(y2in)^==2A-i 
Taking ^ = 1 yields (UMl) 



Ai(e-2/3A) - 5i/3A-i/2Ai'(^-2/3A) ■ 

(371) 
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oo 

e 



C.3. Brownian meander 

We apply again Lemma IC.ll with V given by (|355p ; this time keeping ^ > 
fixed but letting 77 oo] we also assume a; > 0. By monotone convergence, the 
left hand side of (I343|) tends to 

^*E(e-^«X>==(^)'i7(S,(t)); B.,{s) > on [0,i]) dt 

^*]E(e-^«/o ^==(^)^7(i3,(i)) I B,{s) > on [0,<]) P,(to > t)dt, 

(372) 

where Fx{to > t) is the probability that the hitting time tq := inf{t > : 
Bx{t) — 0} is larger than t, for the Brownian motion Bx started at x. For the 
right hand side of l|343p . let us assume for simplicity that f{y) = for all y < x. 
We note that by p56p and (|357p . 4>+{y) does not depend on rj for y > 0, while 
(ji-{y) does. On the half- line [0, cxd), 0_ is by (|342p governed by the equation 
^4>"{x) = (-\/2^a; + X)<j){x); let (/)o and 0i be the solutions to this equation with 

00 (0) - 1, <^'o(0) = 0, (373) 

0i(O) = O, 0'i(O) = l. (374) 

(Note that these do not depend on 77.) Thus, for x > 0, 

,^_(x) = (b-{0)M^) + <^-(0)<^i(x). (375) 

Furthermore, as 77 — > 00, (|357p implies 

0_ (0) Ai(0), (/.'_ (0) ~ -72 771/3 Ai'(O), (376) 

while, by (j358l) . 

«; - -72 77^/3 Ai'(0)Ai(r'/^A). (377) 
As 77 00, wc thus have by p75)) . ((3761) . (pTT)) . 

(?!)-(x) (/)i(a;) _ 01 (a;) 



Ai(e-2/3A) 0+(O)- 



(378) 



Returning to (|343p we thus have, for a; > and assuming f{y) = ioi y < x, 
see also [lOj, Corollary 2.1(ii)], 

-^*E(e-^«/o ^-(^)'*7(B,(t)) I i?,(s) > on [0,t]) P,(ro > t)dt 

'^'^''^ ''^ My)f{y)dy. (379) 



</'+(o) 

Now divide by x and let a; \ 0. The distribution of Bx conditioned on 
Bx{s) > on [0, t] converges to the distribution of the Brownian meander i3me 
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on [0,t]. Furthermore, by the weU-known formula for the distribution of the 
hitting time, with $ the normal distribution function, 

ir,(..>,)4(2*(^)-l)^-J=. (380) 

which also gives the bound 

-P.(ro>t) = i(2$(4)-l) <-|=, (381) 

while the expectation inside the integral on the left hand side of p79p is bounded 
by sup I/I, so dominated convergence applies to the left hand side of (|379p . For 
the right hand side we have (j)i{x)/x <j)'i{0) = 1- Consequently, (|379p yields 



e 







2 



</'+(o) 7o 



4>+{v)f{y)dy (382) 



for every / that is bounded and continuous on (0,cx)) and vanishes on some 
interval (0, <5), and thus by a limit argument (for example a monotone class 
argument) for every bounded / on (0, cx)). Taking f{x) = l/\/2 we find, using 

(|356p and the Brownian scaling Bmc{ts) = t^^'^B^ds), < s < 1, 

Jo VTTi 0+(O)7o .ogoN 

,_,,3 AI(r^/^A) 
^ Ai(e-2/3A) ■ 

Taking £, = 1 yields (fT93l) and A = 1 yields §W^. 

Note also that, since Bmc{t) has the density function {y /t)e^y^ ?/ > 0, it 
follows from p82p that for every y > 0, 



e 







_0+(y)_Ai(^/2el/32; + ^^/3A) 



:dt 



, (384) 



0+(O) Ai(e-2/3A) 
which is a relation for the areas under the Bessel bridges (Bmo | B^c{t) — y). 

C.4- Brownian excursion 

Formula (|384p holds for every A > 0, and by monotone convergence as A ^ 
for A = too. Taking the difference, we obtain 
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e " ' , dt 

V27rt 



Ai(0) 



Ai(e-2/3A) 



(385) 



Let ^ = 1/V2, divide by y and let ?/ — > 0, noting that as y ^ 0, (-Bmo | 

Bmc{t) = y) — ^ i?cx , a Brownian excursion on [0, t] . We then find by dominated 
convergence on the left hand side, and simple differentiation with respect to y 
on the right hand side, the formula ([5T|) . which as explained in Section [Ql yields 
dMl) and HMD. 

Alternatively, we may take f = 1 in (|384p , differentiate with respect to A and 
then divide by y and let y 0. This yields 



d_d_ 
dX dy 



dy 



d_ Ai(\/2y + A) 
y=o d\ Ai(A) 



Ai(V2y + A) _ d N/2Ai'(A) 



Ai(A) 



dX Ai(A) 



(386) 



which yields 



C.5. Brownian double meander 

We do not know how to apply the Fcynman-Kac formula directly to i?dm, but 
we have shown in Section [27] how the result (|314p follows from the result (|310[) , 
or equivalently (|383p . for Smc- 
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